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Abstract: We prove that for a Dirac operator, with no resonance at thresholds 
nor eigenvalue at thresholds, the propagator satisfies propagation and dispersive 
estimates. 

When this linear operator has only two simple eigenvalues sufficiently close to 
each other, we study an associated class of nonlinear Dirac equations which 
have stationary solutions. As an application of our decay estimates, we show 
that these solutions have stable directions which are tangent to the subspaces 
associated with the continuous spectrum of the Dirac operator. This result is the 
analogue, in the Dirac case, of a theorem by Tsai and Yau about the Schrodinger 
equation. To our knowledge, the present work is the first mathematical study of 
the stability problem for a nonlinear Dirac equation. 


Introduction 

We study the stability of stationary solutions of a time-dependent nonlinear 
Dirac equation. 

Usually, a localized stationary solution of a given time-dependent equation rep¬ 
resents the bound state of a particle. Like Ranada |Ran| , we call it a particle like 
solutions (PLS). In the literature, the term soliton is also found instead of PLS, 
but this additionally means that the particle keeps its form after a collision. 
Many works have been devoted to the proof of the existence of such solutions 
for a large variety of equations. Although their stability is a crucial problem 
(in particular in numerical computation or experiment), a smaller attention has 
been deserved to this issue. 

There are different definitions of stability. The first one is commonly called the 
orbital stability. It means that the orbit of the perturbation of a PLS stays close 
to the PLS or a manifold of PLS but does not necessarily converge. A stronger 
notion is asymptotic stability , which means that the perturbation of the PLS 
relaxes asymptotically towards a PLS which is not far from the perturbed PLS. 
In fact in many conservative problems asymptotic stability does not hold. But 
one has asymptotic stability for a restricted class of perturbations, forming the 
so-called stable manifold. 

In this paper, we deal with the problem of stability of small PLS of the following 
nonlinear Dirac equation: 


i dti> = ( D m + V)i/j + VU(i/>) 


(NLDE) 
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where VF is the gradient of F : C 4 i—> K for the standard scalar product of K 8 . 
Here, D m is the usual Dirac operator !Tha92j acting on F 2 (R 3 ,C 4 ) 

3 

D m = a • (-iV) + m/3 = -1 ^ a k d k + m(3 

fc=1 

where m £ K!j_ and a = (aq, tt 2 , a 3 ) and (3 are C 4 hermitian matrices satisfying 
the following properties: 


Here we choose 


where <j\ = 


"E Oi k Oii — 2 5i k 1(£4 , %,k £ { 1,2,3}, 


0 1 
1 0 


b (3ai = 

°c 4 , 

i 

e {1,2,3}, 

1 (C 4 - 




0 CTi\ 

a x 0 J 

and 

p = 

(Ic* 0 
^ o -J c 2 

and 

O' 2 = | 

( 0 -1 
V 1 0 

I and 



In HLDEt , V is the 
such that 


external potential field and F : (C 4 i—> 
V(0, z)elx C 4 , F{e ie z) = F{z). 


K is a nonlinearity 


Some additional assumptions on F and V will be made in the sequel. Stationary 
solutions (PLS) of (INLDEII take the form ip(t, x) = e~ lEt (f>(x) where (j) satisfies 


E<f>= {Dm+ V)(/) + \7F ((/)). (PLSE) 

We prove the existence of a manifold of small solutions to IIPLSEII , interpreted as 
particle like solutions to (INLDE|l . Then we construct a stable manifold around 
this manifold. At the origin, it is tangent to the sum of the eigenspace associated 
with the first eigenvalue and the continuous spectral subspace of D m + V. This is 
the analogue in the Dirac case of lTY02dl Theorem 1.1, non-resonant case]. The 
interpretation is that radiations (described by the continuous spectrum) do not 
destabilize too much the PLS manifold. To prove stabilization towards the PLS 
manifold, we shall need linear decay estimates associated with the continuous 
spectral subspace of D m + V. 

To our knowledge, this is the first stability result on a nonlinear Dirac equation. 


The problem of stability has been extensively studied for Schrodinger and Klein- 
Gordon equations. The methods used to treat these cases cannot be easily 
adapted to our problem, due to the fact that the Dirac operator D m is not 
bounded-below, contrarily to — A. The non-negativity of the latter permits to 
use minimization and concentration-compactness methods to prove the existence 
of orbitally stable standing waves, see e.g. Cazenave and Lions ICL82J or more 
recently Cid and Felrner jCFOll . 

In his review on nonlinear Dirac models, Ranada |Ran| writes that physicists 
first claimed that PLS (Particle Like Solutions) of the nonlinear Dirac equa¬ 
tion couldn’t be stable since the second derivative of the energy functional is 
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not positive-definite. Actually, in a very general setting (not related to the Dirac 
case), Shatah and Straus [SS85j and Grillakis, Shatah and Straus IGSS87I proved 
a general orbital stability condition even if the hessian of the energy functional 
is not positive-definite. Their conditions allow only one simple negative eigen¬ 
value (and a kernel of dimension one also) for the second variation. It therefore 
cannot be directly applied to the Dirac case. However, it gave rise to an inter¬ 
esting discussion about t he a pp licati on of this method to the Dirac equation 
in some physical papers jSV86HAS86IBS V871 . Ranada also refers to numerical 
experiments which seem to confirm that some PLS are asymptotically stable in 
the Dirac case. 

In the Schrodinger case, the asymptotic stability has been extensively studied 
during the last decade. A fundamental work is the one of Soffer and Wein¬ 
stein ISW901SW92| . which is devoted to the study of a small nonlinear per¬ 
turbation of a Schrodinger operator having one simple eigenvalue. They proved 
that the perturbed small PLS relaxes to a PLS. Later, Pillet and Wayne (FW97| 
proposed a different proof in the spirit of the central manifold theorem. In all 
these works, asymptotic stability is a direct consequence of propagation or dis¬ 
persive estimates on the Schrodinger operator. In order to be able to use these 
estimates, one has to to consider the initial state (at time t = 0) of the pertur¬ 
bation as localized i.e. in L l or in L 2 weighted spaces with growing weight. To 
avoid such an assumption, Gustafson, Nakanishi and Tsai |GNT04I proposed to 
use Strichartz estimates. 

Generalizations have been considered for instance by Tsai and Yau ITY02al 
ITY02cirTY02dllTY02bllTsa03| . who treated the case of a Schrodinger operator 
having two simple eigenvalues. An interesting phenomenon appeared: if the two 
eigenvalues are sufficiently distant one from the other, then after linearization 
around the excited state, one obtains a resonance. Tsai and Yau showed that 
if there is no resonance, the manifold of ground state has stable directions. 
In the resonant case, the manifold of ground states is asymptotically stable, 
whereas the manifold of excited states has stable and unstable directions (in 
case of instability, under some conditions, one has relaxation to the ground 
state). For a similar result, see also |SWf)4IISW05| . Notice that earlier Soffer and 
Weinstein [SW99| studied a similar resonance phenomenon in the case of the 
Klein-Gordon equation with a simple eigenvalue; they showed that it induced 
“ metas tability”. Another problem has been studied by Cuccagna (GucOllCucOSI 
IGucOfij . He considered the case of big PLS, when the linearized operator has only 
one eigenvalue and obtained the asymptotical stability of the manifold of ground 
states. Tsai, Yau and Cuccagna also need propagation or dispersive estimates. 
The latter is proved by generalizing the work of Yajima |Yaj95[ on wave operator. 

Interesting development are also given by Rodnianski, Schlag and Soffer |RSS05al 
who proved asymptotic stability of an arbitrary number of weakly interacting 
big PLS. Schlag ISch04j and Krieger and Schlag iKSOfij proved the existence of 
stable direction for unstable big PLS. We point out that some of the works of 
Schlag lKSn4BGS04lRSSn5bl or Soffer lHSS99ll.)SS9lllRSS05bl are dedicated to 
prove dispersive estimates. 

We also would like to mention the works of Buslaev and Perel’mann BEH 
BP92b,BP92c,BP92a , Buslaev and Sulem BS03.BS02 or Weder IWedOOj . in 
the one dimensional Schrodinger case. 
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Here, we study a nonlinear Dirac equation as a perturbation of a linear Dirac 
equation with a Dirac operator possessing only two simple eigenvalues sufficiently 
close to each other. Hence, we avoid problems of resonance after linearization 
around a PLS. The paper is organized as follows. 

In section d we define the important objects and state our main results. We 
start with the propagation and dispersive linear estimates which will be crucial 
tools for this study. Then, we consider the nonlinear equation (INLDEI) and state 
the existence of the PLS manifold. Eventually, we present our main theorem in 
which the stable manifold is constructed. 

The section Q is devoted to the proof of the propagation estimate, which uses 
spectral techniques. This is a time decay estimate in weighted L 2 spaces, express¬ 
ing the fact that states associated with the continuous spectrum are not station¬ 
ary. We use Mourre estimate similarly to Hunziker, Sigal and Soffer jHSS99j (for 
a generalization of the method, see e.g. (BdMGS96] h This method cannot be 
used in the neighborhood of the thresholds which needs a specific treatment. In 
particular, problems can occur in the presence of eigenvalues at thresholds or 
resonances at thresholds, and we shall assume in the whole paper that we are not 
in this situation. For the Schrodinger case, a similar problem has been studied 
by Jensen and Kato QEZ9, Jensen and Nenciu IJN011JN04| . Our arguments 
near the thresholds are inspired of these works. For a related study, see the ar¬ 
ticle of Fournais and Skibsted IfSqH dealing with long range perturbations of 
Schrodinger operators. 

In Section [3 we then prove the dispersive estimate, using the propagation esti¬ 
mate established in Section |5J For an interesting survey on dispersive estimates 
for Schrodinger operators, see Schlag [Sch05| . We have not been able to gener¬ 
alize the methods used in the Schrodinger case, in fact it seems that the Dirac 
equation with a potential behaves like a Klein-Gordon equation with a mag¬ 
netic potential. This fact has already been noticed by D’Anconna and Fanelli in 
[DFj , where they proved simultaneously dispersive estimates for a massless Dirac 
equation with a potential and for a wave equation with a magnetic potential. 
Our method is here inspired of the work by Cuccagna and Schirmer Icsml . 

Finally, the last sections are devoted to the proof of our main result concerning 
the stability of the stationary solutions of lINLDFI . We assume that the Dirac 
operator D m + V have only two simple eigenvalues and that it has no eigen¬ 
values at thresholds nor resonances at thresholds. Note that our assumptions 
exclude electric potentials, for which the theorem of Kramers states that the 
eigenvalues are always degenerate, see IPar9011BH92j . In Section 0] this permits 
us to construct a manifold of PLS and then to study the spectrum of the lin¬ 
earized operator. This in turn, in Sectional will allow us to decompose a solution 
of INLDEII in three parts: the PLS part, the dispersive part associated with the 
continuous spectrum and a part corresponding to “excited states”. This last 
part needs a particular treatment since it is not dispersive and hence disturbs 
the relaxation towards the PLS manifold. 


Acknowledgements. I would like to thank Eric Sere for fruitful discussions and advices during 
the preparation of this work. I wish to thank the referee for useful remarks and suggestions. 
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1. Main results 

This section is devoted to the presentation of the model and the statement of 
our main results. 


1.1. Decay estimates for a Dirac operator with potential. Let us first state our 
results concerning the time decay of e ~ lt ( D m+ v ) i n weighted L 2 spaces and Besov 
spaces. This kind of estimates are called respectively propagation and dispersive 
estimates. As mentioned in the introduction, these results will be very important 
tools for the study of our nonlinear time-dependent Dirac equation. 

The following spaces will be needed to state the main result of this subsection. 

Definition 1.1 (Weighted Sobolev space). The weighted Sobolev space is 
defined by 

H*(R 3 ,C 4 ) = {/ G S'(R 3 ), I m^PYfh < °o} 
for a, t G R. We endow it with the norm 

||/|| h ‘ = \\{QY {P) f fh- 


If t = 0, we write L 2 a instead of Hf. 

We have used the usual notations ( u) = \J\ + u 2 , P = —?'V, and Q is the 
operator of multiplication by x in R 3 . For the sake of clarity, let us also recall 
the 

Definition 1.2 (Besov space). For s G R and 1 < p,q < oo, the Besov 
space ? (R 3 , C 4 ) is the space of all f G <S'(R 3 , C 4 ) (dual of the Schwartz space) 
such that 

H/lk,= (E 2 ^II^*/Hp) <+°° 

\ieN J 

with (p G 2?(R ra \ {0}) such that (5(2 _:, £) = 1 for all (GR 3 \ {0}, <Pj{Q = 
t for all j G N* and for all £ G R 3 , and tpo = 1 — Tj- We endow it 

with the norm f G Bp q (R 3 , C 4 ) i—> 11/11s® g - 

In the whole chapter, we shall work within the following 

Assumption 1.1. The potential V : R 3 i—> 64 (C) (self-adjoint 4x4 matrices) 
is a C°° function such that there exists p > 5 with 

Vo G N 3 , 3 C > 0, \/x G R 3 , \d a V\{x) < ——— T . 

(a;)p+l Q l 

Notice that by the Kato-Rellich Theorem , the operator 

H := D m + V 

is essentially self-adjoint on C“(R 3 , C 4 ) and self-adjoint on H 1 ( R 3 , C 4 ). We also 
work with the 
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Assumption 1.2. The operator H presents no resonance at thresholds and no 
eigenvalue at thresholds. 

A resonance is an eigenvector in H^ 2 (M. 3 , C 4 ) \ C 4 ) for some cr £ 

(1/2, p— 1/2) here. Let 

P c(H) = t { — oo,—m]U[-t-m,+oo) (H) (1-1) 

be the projector associated with the continuous spectrum of H and 

Tt c = P c (H)L 2 (M. 3 , C 4 ). (1.2) 


We are now able to state our 


Theorem 1.1 (Propagation for perturbed Dirac dynamics). Assume that 
Assumptions Ell and El hold and let be a > 5/2. Then one has 

\\e-' tH P c (H)\\ B(LltLlir) <C(f>- 3/2 . 


The proof of this result will be given in Section El We notice that it is still true 
if we assume p > 3 in Assnmption ll.l| 

Our next result is the following theorem, proved in Section E 

Theorem 1.2 (Dispersion for perturbed Dirac dynamics). Assume that 
Assumvtions M. 1\ and il. A hold. Then forp £ [1,2], 9 £ [0,1], s—s' > (2+0)(| — 1) 
and q £ [l,oo] there exists a constant C > 0 such that 


||e- ltff P c (iO|| B . <C(A(t))i " 1 

P><2’ p',q 

with - + A? = 1. and 

p p' ’ 

f \t\~ 1+e/2 if\t\ £ (0,1], 
K(t) = l 

{ \t\ 1 6/2 if \t\ £ [1, oo). 


1.2. The stable manifold around the PLS for the nonlinear Dirac equation. We 
now want to study the following nonlinear Dirac equation 

f i dt-ip = Hif + VF(V>) n 

U(0 ,-)=^o. 


with if £ C 1 (/, H 1 ( R 3 , C 4 )) for some open interval / which contains 0 and where 
we recall that H = D m + V. The nonlinearity F : C 4 h> R is a differentiable 
map for the real structure of C 4 and hence the V symbol has to be understood 
for the real structure of C 4 . For the usual hermitian product of C 4 , one has 

DF(v)h = $t(VF(v), h). 


We work within the following 
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Assumption 1.3. The operator H has only two simple eigenvalues Ao < Ai, 
with <j )o and (f >i as associated normalized eigenvectors. Moreover, the non reso¬ 
nant condition 

|Ai - A 0 | < min{|A 0 + ro|, |A 0 - m\} 


holds. 


Assumption 1.4. The function F : C 4 i—> R is in C°°(R S ,R), is a homogeneous 
polynomial of degree 4 (i.e. with D a F(z) = 0 for |ct| = 5 and D@F( 0) = 0 
for \/3\ < A) or satisfies F(z) = 0(|,s| 5 ) as z — > 0. Moreover, it has the gauge 
invariance property: 

F{e w z) = F{z),\/z £ C 4 , V9 £ R. 

We will prove in Theorem roi that some solutions of the equation in are 
global and can be decomposed as the sum of a PLS plus a remainder part which 
is vanishing. Since the PLS part may change during the evolution, we need to 
track it. So we prove that around the origin, PLS form a manifold. We have the 


Proposition 1.1 (PLS manifold). Suppose that Assumptions I /. /I| 1 ,J\ hold. 

Then for any a £ R + , there exists FI a neighborhood of 0 £ C, a C°° map 

h-.Fl {</o} X n H 2 (R 3 , c 4 ) n Ll(R 3 , c 4 ) 

and a C°° map E : FI i—> R such that S(u) = u(j >o + h{u) satisfy for all u £ FI, 
HS{u) + VF(S(u)) = E{u)S{u), (1.4) 

with the following properties 

( h.(e ie u) = e l6 h(u), VO £ R, 

i h(u)=o(\ u n 

| E{u) = E(\u\), 

{E(u) = X 0 + O(\u\ 2 ). 

Proof. This kind of results is now classical and left to the reader. For more 
details, see Subsection FHI 

We are now able to write the main theorem of this paper. Its proof is given in 
Section 6 . To state it we need the space H c defined in n~2i 

Theorem 1.3 (Stable manifold). Suppose that Assumptions I /. hold. 
Let s,s',/3 £ Rlj_ be such that s' > s + 3 > /3 + 6 and a > 5/2. There ex¬ 
ists £o>0, 1?>0, K >0 and a Lipshitz map 

F : B c (0, e) x (h c n B h .j (0, R)^j ^ C 

with F(v,0) = 0, 

|^,0l<A'(|f;| + ||e||^) 2 , 

and such that the following hold. For any initial condition of the form 
ifo = S{v 0 ) +Co + &(vo,£o)(l>i 
with vo £ Bc(0,e) and € TL C D B hs i (0, R), one has 
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(%) there exists a unique global solution ^ of (jl .31) in 

C°° (K,H s '(R 3 ,C 4 )nifl CT (IR 3 ,C 4 )n J B^ i2 (R 3 ,C 4 )) ; 

(ii) there exists (voo ;£oo; Eoo ) £ C x H* nH c x R with 

K'oo - voI < K\M 2 H s', |£ool < k\m 2 h .„ HCoo - £oll/fa' < k\M 2 h ,, 

Ip(t) = e~ l( ' tE( ' Va °^ +E °°^ S(voo) + e~ ltH ^ao + e{t), 

< MQWhl. < ^II&IIhs' 

IK*)II< 2 — Jj^2 H^ollffX- 

Remark 1.1. The proof of these theorem work also if we want to obtain an ex¬ 
pansion of the form 

m = e~'( tE ^ +E ^S( Voo ) + e ~ ltD + e{t) 
with the free Dirac operator. But in this case, we only have the estimates 


such that 

where 


as t —> +oo. 


Coo Co 


H 


< ^IlCollff. 


see 


the remark following Lemma OB 


We notice that the stabilization is “faster” than the propagation and the dis¬ 
persion: it is of order ( t)~ 2 whereas e _ltff Coo is of order (t) _3//2 by Theorems ll.il 
and ll.2l Hence the theorem states the existence of a family of initial states which 
form a manifold tangent at the origin to the sum of the eigenspace of H asso¬ 
ciated with Ao and the subspace associated with the continuous spectrum of H : 
H c - This family of initial states gives rise to solutions of Cl which asymptoti¬ 
cally split in two parts. The first one is a PLS: e - 1 ( t£: (“oa)+.E 00 ) S(uoo) the other is 
a dispersive perturbation: e _lt ' f4 Coo- Hence if one perturbs a PLS in the direction 
of the continuous spectrum then this PLS relaxes to another PLS by emitting a 
dispersive wave. 

This phenomenon is due to the propagation and the dispersion properties of the 
subspace associated with the continuous spectrum of H. We don’t think that 
such a phenomenon could take place for perturbations in the direction of the 
excited states <j>\. Indeed, on this subspace, the dynamic seems to be conservative. 
The fact that we use propagation and dispersive estimates restricts the family 
of perturbations to regular and localized ones. 

We now turn to the proof of our results. 
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2. Proof of Theorem ll.lt propagation estimates 

Here we prove the propagation estimates of Theorem 11.11 The method used 
by Jensen and Kato to prove this kind of estimates for Schrodinger 

operator works only for initial states which are spectrally localized near the 
thresholds ±m. They used the spectral density as the Fourier transform of the 
propagator. But the Dirac resolvent 

R v ( A±ie) = (H- A Tie )” 1 

does not decay in B(L^, L 2 _ tT ) as |A| — * +oo for any er > 0, see ES m Ml- So we 
cannot use its Fourier transform. To our knowledge, this method is the only one 
that permits to treat the problem of propagation for energies near thresholds. 
Hence with this method, we only prove (in the section ITU the 

Proposition 2.1 (Propagation near thresholds). Suppose that Assumptions 
to and roi hold and let \ € CjJ°(R 3 ,C 4 ) be such that its support is in a suffi¬ 
ciently small neighborhood of [—m; m]. Then one has for a > 5/2 

\\ e ~ ltH Pc (H) X (H) II B(Ll,L*_,) < C (t)~ 3/2 . 

We recall that P C (H) is defined by 11 .1 1) . 

We also need to treat the propagation estimates for initial state whose spectrum 
does not contain any threshold. We cannot use the spectral density. So we work 
directly with the propagator. This is exactly the method used by Hunziker, Sigal 
and Soffer in 1HSS99I . B ut in our case, their result needs some adaptation. Hence 
we need to generalize 1HSS99I Theorem 1.1] to the case of unbounded energy. In 
Section EM we prove the 

Proposition 2.2 (Propagation far from thresholds). Suppose that Assump- 
t.ian l I. H hnlds. Then for any \ G C°°(R 3 , C 4 ) bounded with support in R\(—m; m) 
and for any a > 0, there is C > 0 such that 

\\e~ ltH x{H)\\ B{Ll , L *_ a) <C(t)-°. 

The proof of Theorem o is then a consequence of the above propositions 
Proof (Proof of Theorem E W - 

We choose xo G C°°(R 3 , C 4 ) satisfying the assumptions of ProDQsition l 2 .il Xoo € 
C°°(R 3 ,C 4 ) satisfying assumptions of Proposition 12.21 such that xo + Xoo = 1- 
Hence the continuous spectrum of H is divided in two parts. We obtain the 
inequality 

l|e- ltH Pc (H) < \\e~' tH Xo (H) P c (H) \\ B{L *, LU) 

+ ||e ltH X oo (H)\\ B{LliL i a y 

Hence from Proposition rm and EM we deduce Theorem o 
It therefore remains to prove Propositions EM and EM 


2.1. Step 1 : Propagation near thresholds. 
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2.1.1. Proof of Proposition \2.1\ We now prove Proposition rm Let x be in 
Cg°(R 3 , C 4 ), then the operator e~ ltH P c ( H ) % ( H ) as a function of t is the Fourier 
transform with respect to A of 

A ^ QB+( A)l ( — oo,—m]u[m,oo)(A),\(A), 

where 

Ry( A) = lim i?y(A±ie), (2.1) 

E—>0+ 

we will prove in Section lT^l that the limit exists in B(LL 2 _ a ). So Proposition l2.il 
is a consequence of the 

Proposition 2.3. Suppose that Assumvtions il . il a,nd \T~B hold. Then for A > m 
close enough to m, one has 

Ry( A) = lim AV(A ± ie) 


exists in B(R a 1 ^ 2 , for a > 3/2. It is C l if a > 1/2 + 1 and 0 < l < 2 with 

A) = 0{V\- rn /2 l ), (2.2) 

as A —» to + . 

T/ie same holds for A < —m */m is replaced by —m. 

We prove it in Section 12. 1.21 The idea is then to apply to 

A t—> 9i?y(A)l(_ 00 j _ m ] U [ TO> 00 )(A)x(A). (2-3) 

with k = 1 and 9 = 1 / 2 , the following 

Lemma 2.1 (Lemma 10.2 of . TK79 ). Suppose F( A) = 0 for A > a > 

0, F'( fe + 1 ) £ L 1 {{5 1 +oo[) for any 5 > 0 and an integer k > 0 and that F^ k+1 ^ (A) = 
0(X e ~ 2 ) near 0 for some 9 £ (0,1). Assume further that F^l (0) = 0 for j < k—1, 
then 

F(t) = 0{t- k ~ B ). 

The symbol O may be replaced by o throughout. 

We refer to rnrroi for the proof of Lemma rm In fact to apply this lemma 
to (12.311 . one should split this function in two parts, one supported in R + and 
the other in M~. Then one translates the first one by — to and applies the lemma. 
To deal with the other part, one works exactly in the same way after a symmetry 
with respect to the origin. To end the proof of Proposition ^. II it remains to prove 
Proposition 12.31 This the goal of the next section. 
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2.1.2. Behavior near thresholds of the Dirac resolvent: proof of Pronosition WI\ 
In this section, our aim is to prove Proposition roi First of all, we notice that 
if the limits ED exist then we have 

Ry(\y = R+(\), 

and since 

<*5 (An + V - z)~ 1 Q'5 = -(Dm + a 5 Va 5 + z)" 1 , 

3 

for c ^5 = cti/3, one obtains 

i= 1 

a 5 R±(\)- l a 5 = -Rl Va 5 (-V- 

So we only need to study the behavior of Ry{ A) near +m. Let us introduce 

C++ = {z £ C, Qz > 0, 5iz > 0} 

then the behavior for the free case (V = 0) is given by the 

Proposition 2.4 (Dirac’s resolvent expansion). Let be s, s' > 1/2 with s + 
s' > 2 and t £ R. Then Rq(z) £ H[_ s ,) is uniformly continuous in C++ 

and so it can be continuously extended to C++. 

Moreover, the formal series z £ C++, 


Ro{z) = z 2 - m 2 yD m Gj + Y, Z ( 1 Vz 1 ~ m 2 yGj 

3=0 j=0 


with 3(\/ z 2 — m 2 ) > 0 , is an asymptotic expansion for z —> m in the following 
sense: 

Let k £ N, if Rq{z) is approximated by the corespondent finite series up to j = 
k, the remainder is o(|z — m\ k / 2 ), as z —> m, in the norm of B (Hl^ 1 , Hf s ,) 
with s, s' > k + 1/2 (and s + s' >2 if k = 0 / and t £ R. 

In the same sense, this identity can be differentiated in z any number of times. 
More precisely, for l £ N* the I th derivative in z of the said finite series is equal to 
■jp-R(z) up to an error o{\z — m\ k / 2 ~ l ), as z —> m, in the norm ofB (Hf^ 1 , Hf s A 
with s, s' > k + l + 1/2 and t £ R. 

Proof. It is an adaptation of lemmas of |JK79J. We rewrite nTTrol Lemma 
2-1], IJK791 Lemma 2.2] and . .1K79I. Lemma 2.3] in the Dirac case with help 
of the identity 

{D m - z)~ 1 {D m + z ) _1 = {-A + m 2 - z 2 ) -1 , 


or in (C 2 ) 2 


{Dm - z) 1 = 


z + m 


a ■ V 


—A — z 2 + m 2 — A — z 2 + to 2 
a • V z — m 


\ -A - z 2 + to 2 —A — z 2 + m 2 


where a are the two dimensional Pauli matrices. 
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To obtain the behavior of the Dirac resolvent in the general case, we would like 
to use the formula 

R v {z) = M(z)~ 1 Ro(z ) (2.4) 

with 

M(z) = (l + R 0 (z)V). 

To give a meaning to Identity (IP) . we have to prove that M(z) is invertible 
in B(H]J 2 ) for a > 1/2 with er+1/2 < p , where p is introduced in assumption ll.il 
We will also give the asymptotic behavior of Ry(z) and some of its derivatives 
as A —> m + . By means of Proposition 12.41 one has 

2 e C77 1 ► M(z) e B (H 

is uniformly continuous for 1/2 < a and 2 < a + o' < p and some a' > 1/2. We 
also have 

M (z) = M (m) + A(z), 

with A(z) uniformly continuous in B ^H 1 / 2 ,, H 1 / 2 ^ near m in C++ and tending 

to 0 as A —> m for 1/2 < a and 2 < a + a' < p and some a' > 1/2. We now 
prove the 

Lemma 2.2 (Threshold’s eigenvector and resonance). Suppose that As¬ 
sumption E3 holds. Let A4(s) be the kernel of M(m) in H l J_ 2 and IC(s) the 
1/2 

kernel of (H — m) in H_ s . Then A l(s) and /C(s) are finite dimensional and do 
not depend on s £ (1/2, p— 1 / 2 ). So we write Ai and K. and we have 

M = JC 


Proof. See also rnpi Lemma 3.1]. 

Let u £ IC(s), then ( D m + V — m)u = 0 and u £ so Vu £ H~^ s 2 

and since p — s > 1/2, s > 1/2, and s + p — s > 2, we obtain, by Proposi¬ 
tion [^1 ( D m — m) 1 ( D m — m)u = ( D m — m) 1 Vu £ H]/ 2 . For any <f> £ C“, 

(, (Dm - m.y 1 (D m - m)u) = ((D m - m) (D m - my 1 0, u) = (</>, u), 

we obtain (D m — m)(D m — my 1 Vu = Vu and (D m — m)(u + (D m — m) 1 Vu) = 
0. Since D m — m has no kernel in Hj s , because there’s no harmonic function 
in L 2 _ s , we obtain u + (D m — m ) -1 Vu = 0. Hence, we have K(s) C A i(s). 
Conversely, I + (D m — m ) _1 V defines a Fredholm operator of B(H]/ 2 ). If u £ 
M(s) then u € H]/ 2 and (D m —m)~ 1 Vu £ H 1 / 2 . So we write 0 = (D m — m) (n+ 
(D m — my 1 Vu) = (D m — m + V)u and we obtain M(s) C M(s). 

Now we introduce I + V (D m — m ) -1 £ B^f 1 ^ 2 ), and its kernel Af(s) which 
is finite dimensional is a Fredholm operator. We have that A f(s) is decreasing 
with s and A4 (s) is increasing. Since, by duality, dimA4(s) = dimA/"(s), we 
deduce that J\f(s) and K.(s) = M (s) do not depend on s. 

We are now able to conclude the proof of Proposition ^. 31 
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Proof (Proof of Proposition I Assumption I I .21 gives /C = 0 and so with 
Lemma IT2l one obtains M = 0. Hence M(m) is invertible since it is a Fred¬ 
holm operator. We use Von Neumann series to obtain that M(z) is invertible 

in B ^ ft]/ 2 j for a > 1/2, 2 < a + a' < p and some a' > 1/2 and 

M(z)- 1 = M (to ) -1 J2(A(z)M(m)- 1 ) j . 
j> o 

So for A > m close enough to m, 

M+( A)" 1 = lim MiX + ie)- 1 
£—>■ 0 + 

exists in B (h]/ 2 ^ for cr > 1/2 with 2 < cr-|-cr' < p and some a' > 1/2. We obtain 

that lim Ry (A + ie) exists in B (h^^ 2 , H^ 2 \ for cr > 1/2 and cr > cr" >1/2 
£—>0+ V / 

with o + cr” > 2,2 < a + cr' < p and some o’ > 1/2. 

Using Proposition ^. 41 we prove that if 1/2 + fc < cr, and a' + 1/2 + k < p then 

fjjk M+ ( A) = 0 (VA-m 1/2 fc ) 

in B [h_ 1 J 2 , H 1 / 2 ^ for fc £ N* as A —> m + . Since we have 

Tx fix )"‘ = (J;*™) f(a >" 

for matrix valued differentiable function U with invertible values, we obtain 
for fc G N* the estimate 

^M+(A) _1 = 0(VA - m 1/2 fe ), 

in £> with 1/2 + fc < a and cr + l/2 + fc<pasA—► m + . So by Leibniz 

formula, we also have for fc G N* 

f[Xk R vW = 0(V\-m 1/2 k ), 

in B {pt a }! 2 , with 1/2 + j < cr, 1/2 + fc — j < cr' and 1/2 -|- fc — j + a < p 

for all j G {0,..., fc} as A —> m + . For the case fc = 0, we have the formula 

Rv(z) = Rq(z) (1 + VRo(z ))~ 1 . 

Since R + (m) = R~(m), this leads to 

SsRy(m) = 0, 

and so 

3tR±( A) = 0(V A“^ 1/2 ), 

as A —* m + in B (^Ha 1 ^ 2 , H^ 2 ^j with 3/2 < a and cr + 3/2 < p. Hence (12.21) is 
proved. 
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2.2. Step 2: Propagation far from thresholds. In this section, we prove Proposi¬ 
tion P We prove the propositions for t > 0. Then using (e ltH y = e ltH , the 
result easily follows for t < 0. 


2.2.1. Proof of Proposition 12.21 Let us introduce 

a = ±{d- 1 p-q + q-pd- 1 }. 

Lemma 3.1] gives that A is an essentially self-adjoint operator and the 
domain of its closure contains the domain of ( Q). Proposition 12.21 is then a 
consequence of the 

Theorem 2.1 (Minimal escape velocity). Suppose that Assum,ntion \l. 1\ holds. 

Then for any x £ Cff bounded with support in (—oo, — m ) U (to, +oo), there ex¬ 
ists 9 > 0 such that for any l £ R, for any v £ (0,0), and any a £ R one 
has 

Vt > 0, \\t A -a-vt<Oe~ ltH X( H )^A-a>o\\ < Ct~\ 
where C do not depend a and t. 

The proof will be given in Section 12.2.21 Let us now show that Theorem rm 
implies Proposition ^. 21 

Proof (Proof of Proposition QO|l. We notice that for c > 0 
{A)~ a = (A)~ a t±A>ct + 0(t~ a ), 
when t > 0, this leads to 

(■ A)- a e~ itH x(H)(A >"“ = (d)-“l 4<1 M,e-“ X (/l)l A> «(i)-“ + 0(r”), 

2 — 2 

So if we choose a = — ^ and v = 9 — | in Theorem l2.ll we obtain 
\\{A)- a e- itH x(H)(A)- a \\ < Ct~ min(a ’ l) . 


Then we prove that ( A) a (Q) “ is bounded for any positive a. It is quite immedi¬ 
ate for integer a using multi-commutator expansion iHSOOl Identity (B.24)]. To 
prove it for any positive real, we use ; SS'981 Identity (1.2)]. This identity states 
that for a self adjoint with B > 1 and a positive real 3, we have on domain 
of BM+ 1 


B , =S n 

7T J o B + W 


where {/3} =/?—[/?] and [/?] is the integer part. With this formula for B = ( A) 2k 
for any k £ N, we prove for any (3 £] 0, 1[ that 


{A) 2k0 < C{Q) 2k0 . 


This ends the proof of Proposition 12.21 
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2.2.2. Proof of Theorem, 12.11 Our proof of Theorem rm is an adaptation of the 
one of JHHMS], we make some modifications. 

For any self-adjoint operator B with domain D(B), we write AcIa{B) for the 
operator [A, B] with domain D(A) D D(B) dense in D(B ), defined by 

Vrt, v £ D(A) (~l D(B), (i[A, B]u, v) = i ((Bu, Av) — ( Au , Bv)). 

First of all, we have 

Lemma 2.3. Suppose that A ssv.m.nt.inn i 1 . 1 i holds. Then Ad\(H) is bounded and 
can be written as a finite sum of terms of the form 

f(P)g(Q)h(P) 

where f and h are rational fractions with coefficients in Al^C) of degree at 
most 0 with no poles, and g is a function that satisfies Assumption \l.ll 

Proof. The proof is a simple calculation based on the fact that Adp-(f(Q)) = 

- 1 (Vjf)(Q). 

We can state the 

Lemma 2.4 (Mourre estimate). IfV satisfies Assumvtion ll . 11 then for any 
9 £ (0,1) there exists v > 0, one has 

and for any A £ (— 00 , —m) U (+m, + 00 ) for any 6 > 0 there exists e > 0 such 
that 

^]1|H-A|<e > ( 9 9 = _ l|ff-A|<e- 

VvA^-l -m z J 

Proof. This is a consequence of 1 [D m ,A] = zrsWrP anc ^ a com P ac t 

operator in B(L 2 (M. 3 , C 4 )). 

We now adapt [HSS99I Theorem 1.1] to the case of unbounded energy since 
here the multi-commutators Ad k A {H ) are bounded operators. We need the We 
introduce the 

Definition 2.1. We call generalized indicator function of R~ a function of the 
form 

u(x) 

x 1—> e -5- 1 R - (x) 

with u £ Cq° (R) supported in [—77,77] (for some 77 > 0 ), nonnegative, and such 
that 77(0) = 1 . 

Note that our generalized indicator function of R _ are of infinite order in 
the sense of IH5'S99l , Section 2]. Using the commutators expansion presented 
in IHS00I Section B] and the Mourre estimate of Lemma \2 .41 we have the 
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Lemma 2.5. Let f be a generalized indicator function of Mr and g £ C°°(R) 
with support sufficiently far from thresholds ±m or support sufficiently small in 
(—oo, — to) U (+m,+oo). Let be A s = s _1 {T —a} and 0 < e < 1. Then for 
any n £ N and 5 > 0 there exists a C > 0 independent of a £ R such that 
for s > 1 


g(H)i[H,f(A s )]g(H) < S - l 0g(H)f (A s )g(H) 

+ Cs~ 1 ~ e g(H)f 1 ~ s (A s )g(H) + g 2 (H). 

Proof. See [BSS9Q1 Lemma 2.1], in our case we don’t need to replace H by 
b(H)H with b £ Cf°. Indeed, our commutators Ad\(H) are bounded by means 
of Lemma Q Then we replace the notion of function of order p by the one 
of generalized indicator function. Finally, we use the fact that a generalized 
indicator function / satisfies 


Wk£ N, V<5 £ (0,1), 30 0, 


fW <c |/| 


1-5 


We are now able to give the 

Proof (Proof of Theorem. 1 2. II) . We write y as a finite sum of function gj £ C°° (R) 
with support sufficiently far from thresholds ±m or support sufficiently small in 
(—oo, — m) U (+m, Too). If we prove the theorem for gj instead of x the theorem 
follows by summing each estimates for gj since the sum is finite. In the rest of 
the proof, we will not write the index j of g. 

We notice that ifO<u<0 — g and if F is a positive non increasing C^-function 
which equals 0 on R + and 1 on (—oo, —rj), we have 

l(A-o-vs)<0 < F ^—-- 9 j ■ 


Now suppose F is a generalized indicator function of R . We consider 

F(s~ 1 {A-a-et}) 

and study the time evolution of the observable g(F[)f(A ts )g(H), where / = F 2 , 
with respect to e~ ltH tA-a>o- That is to say we study 

{e~' tH lU-ax,?/’, g(H)f (A ts )g(H)e~ ltH l A _ a > 0 V’) • 

We work exactly as in the proof of [HSS991 Theorem 1.1]. Hence using Lemma ITT)! 
we obtain for 0 < t < s and s > 1 


(e- ltH t A _ a> 0 fj,g(H)f(A ts )g(H)e- ltH t A - a >oi() < C S -^ n ~ 2 -^\\fj \\ 2 

PCs- 1 -* [ (e- I ^l A _ Q>0 ^,ff(i?)/(A ts )ff(i?)e- ltff l y 4_ a >oi/>) 1 - 5 ||^|| 25 . 

Jo 

Then using the Gronwall’s lemma (see |ABdMG96ll Lemma 7.A.1]), we obtain 
{e~ ltH 1 A-a>o^, g(H )f{A ts )g{H ) e ~ ltDm 1 A - a>0 if) 

< {Cs- s{2n - 2 - £ )\\i(\\ 25 + 6Cs- e \\i/j\\ 2s y /5 , 
so if we choose a small S and a big n, the proof is done if we choose s = max {1, t}. 
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3 . Proof of Theorem 11.21 dispersive estimates 

Dispersive estimates for Schrodinger operators with electric potentials take place 
in Lebesgue spaces. This fact permits to use simple perturbation methods (like 
Duhamel’s formula) to prove the decay estimates for perturbed Schrodinger 
equations. Unfortunately, we have only been able to prove dispersive estimates 
for Dirac operators in Besov spaces, so it was not possible for us to use Duhamel’s 
formula or other perturbation method used for Schrodinger operators. 

We notice that in the case of a Dirac operator with scalar potentials (matrix 
valued functions colinear with /?), the square of the Dirac equation gives four 
coupled Klein-Gordon equations with an electrostatic potential. This permits to 
use results on the Klein-Gordon equation. For example, Yajima |Ya^5| proved 
dispersive estimates for the Klein-Gordon equation by using wave operators as¬ 
sociated with Schrodinger operators including an electrostatic potential. But in 
the general case, by taking the square of a Dirac operator with a potential, we 
obtain also a magnetic potential. Hence the method used by Yajima does not 
work in our case. 

To our knowledge the only one study of the dispersive estimates associated 
with the Dirac equation, is the the work of D’Anconna and Fanelli IPFI for the 
massless case. For non zero mass we have not been able to found any reference. 
Even for the free case for which dispersive estimates can be deduced from those 
of Klein-Gordon equation. Here, to give a sketch of the proof for the general 
case, we first prove the free case estimates (see Section E3, using estimates 
on oscillatory integrals of Section IXTI In Section El.31 following Cuccagna and 
Schirmer IcsoTl . we introduce the distorted plane waves. This permits us to 
tackle the proof of the general case in Section Hi. 3. 21 


3 . 1 . Estimates on some oscillatory integrals. Here, we state some stationary 
phase type results which will be useful for the rest of the proof. We denote 
by S 2 the unit sphere of R 3 . 

Lemma 3 . 1 . Let be f £ C 1 (S' 2 ) and for any v £ S 2 and any k £ R define 



Then we have 



(3.1) 


where C does not depend on f , k or v. 

If f is in C 2 (S 2 ) with f(v) = f{—v) = 0, we have 



(3.2) 


where C does not depend on f , k or v. 
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If f is in C 2 (S 2 ) and vanishes in a neighborhood of v and —v, we have 


\Jy{k)\ < 


C 

( fc ) 2 


|V 2 “/(u 
M<lJS 2 \u-v\'° 


S|a|<l ft 


■duj J2 


dist(supp(/), v) 


+ 


f |v 2a /M ' 
Js 2 duj 


M<1 JS 2 |aj+t?| 


dist(supp(/), -v) 


(3.3) 


where C does not depend on f, k or v. 

Proof. We can suppose v = (0,0,1) since estimate (13.111 . (13.211 and 13.31 1 are 
invariant under the action of rotations. We have 

/• 27T /»7T 

J v (k)= / 0) s in(</>) d<t>d6, 

Jo Jo 

then we make an integration by parts in </> 


J v (k ) = - 


O ik{l—cos((f))} 


m<t>) 


dd 


J o 

i 


f‘2ir /*7r 


D l/c{l —COS' 


W»fy/(M) d4dS, 


If we suppose that / vanishes in a neighborhood of v or —v, then we use that 
for any qY 

\f{0A')\ < f|*/(MI 

Jo 

to obtain en> in this case. Otherwise with help of a smooth cut-off, we split the 
integral in two parts, each one has a support far from v or —v. Repeating the 
previous proof for each part, we prove the estimate in the general case. 

If moreover we have f(v) = f(—v) = 0 then we have for any a > 0 


p2ir not 

J v (k) = -j J g^U-cosW)} d^f{9,<ji)d4dB 

/* 27 T /» 7 T 

+ t / e lk ^- cos ^d^f(e,ci>)d(pde 

k Jo Jn-Ot 

n2n n'K—Ot 

+ J / gifco-cos m d<k f(e,ct>)d(i>de. 

k Jo Ja 

We use an integration by parts to obtain for the second term of the right hand 
side 



e iMl-cos 


0 \k{l— cos(0) } f (^ 

sin ((/>) 


i f V “ lfc {i-cos(0)} f ^/(M) _ cos(^)d < j ) f(8,4>) 

k J a [ sin(0) sin(c^) 2 


d<p, 
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for the other terms of the right hand side direct estimations give us 

\Jv(k)\ < 777 / sup|<V(M)l d6 
\k\ Jo <t> 

S'1 /*27T [ /»27T /*7T 

+ _ Tid2 / \ su Pl«W(M)l d6+ / / |d$/(0,0)| dcfdd 

a \ K \ Jo 4> Jo Jo 

choosing a = \J\k\ and working like in the proof of the estimate ED , we 
obtain estimate (IQ) . The reader recognized the proof of the well known Van 
der Corput Lemma with modification in order to give precise estimates. 

For the estimate iTHli . we first split the integral J v (k) in two hemispheres with 
respect to the pole v and we choose a = dist(supp(/), v) or a = dist(supp(/), —v). 

We obtain first the 

Proposition 3.1. Let h £ C(R) and g £ C 2 (R 3 ) be such that the integrals ap¬ 
pearing in the following estimate are finite. Then defining 

I{k,u) = [ 

Jr 3 


for any u £ R 3 and any k £ R, we have 

il«M l ? °3«>l dt [ | £ |M-i_!^»(«)I 

J R 3 


\I{k,u)\ < r^7 max i [ |£|' 
\ku\ |a|<l Jr 3 


m m 


M 


U I £ 

hi + Id 


hi 


where C does not depend on h, g, k or u. 

If moreover g vanishes in a cone of axis D = Span(u), we have 


df 


(3-4) 


\I(k,u)\ < 


C 


\ku\ dist(supp(g) D S' 2 , D fl S 2 ) 

2 IV 2 “3(0' 


x max 
H<i 



2|q| — 2 _ 




M 


df, [ |^| 2|Q| - 2 

J R3 


— 2 V 2 “ ff (0 


i«i ici 

where C does not depend on h, g, k or u. 
Proof. We write 


hi ^ Id 


hi 


df 


I(k, u) = [ e i*=WI«!)-*■“}<,(£) d£= [ e lk ^-^Ju^pk\u\)p 2 dp, 
J R3 J R+ 1 1 


where J v , P (k) = f s2 e 1 *^ 1 v '^g(puj) dw and we apply Lemma IXTl 
We introduce a first useful variant with the 
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Proposition 3.2. Let g G C 1+fc (R 3 ) be such that the integrals appearing in the 
following estimate are finite. We introduce 


F(x ) = f e 1 

J R 3 




} g(0df 


for any i£l 3 . Then for all a £ N 3 such that |a| < k we have 


all afN 3 such that |a| < k, we have 


|V^<?(£)| 

-dA. 

(3.5) 

x £ 

l*| |f| 

J 



Then for 


\S7 a F(x)\ < 


C 


;r|l“l +2 dist(supp(g) n S 2 , D+ n S 2 ) 


x max 

\0\<2+\a\ I J R 3 


| ? | 1 / 31 - 2 ^ 5 ( 0 ^ ^ . ( 3 . 6 ) 


R HI 


Proof. The critical points correspond to the the semi axis spanned by x. We treat 

the part of the integrals which is far from critical points by using an integration 

£ — .rc 

by parts with help of the operator L = ^ ^ ^ • Vt. Let be <j>(x, £) = {|£| |.t| — 

RlTeT _ FTl 

£ • a;}, we have 

F(x) = (Le^’-\g) = (e'^\L*g), 

with 

2 

L* = -L --- 


Ml£l 


ill R 


This gives the bound 


n^ax{/ |£|I /3 I _1 |V ,3 5(£)| df \ , 

1*1 \P \< 1 l Jr 3 J 


or after an iteration 

02 ma/ |£|^ _1 |V^g(£)| df \ , 

1*1 ldl<2 (J R3 J 

we obtain Estimate 13.51) for a = 0. The method to treat the other part of the 
integral is exactly the one we used in the proof of Proposition 13.II 
For higher order derivatives, we have 




1*1 


and so 

(V x e^),g) = -±-(e'^’-\X7\Q\g). 

1*1 

the result is then obtained by applying this trick |a| times and then repeating 
our proof for the case a = 0, we obtain Estimates ifTol) and ESJ for X7 a F(x). 
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And finally, we need the 

Proposition 3.3. Let be g £ C 2 

R 3 , k £ R and 


we have 


C 


\k\ 


3/2 


with compact support. Then for any u £ 

I(k,u)= [ e ^{V? 2 +™ 2 -^} 

J R3 


(0'“ hl v a fl(0 


max 

max < 

I 


M<2 l 

J R3 


inf 


: if ]L m 


l—n—1 


0f £ ,«O 


«_ L 

I«1 151 


df 


(3.7) 


\ xesu PP(s) 1 'Jx 2 + rri^ 


Proof. We can suppose u = (0,0, |it|) since estimate (13.71) is invariant under the 
action of rotations. 

The oscillatory integral I(k,u) is bounded and critical points of the phase 
of I(k,u) are supported by the semi axis spanned by u. With help of a smooth 
cut-off function \ > we split the integral in two parts J(fc, u) = Ii(k, u) + l 2 (k, u), 
where I\{k,u) is supported in a half cone around u. We then use multiple inte¬ 
grations by parts with help of the operator 


L = 




— U 


1*1 


yj £ 2 +m 2 


■V«. 


Since (1 — x) 9 £ C 2 (R 3 ) has support far from critical points and since for A(£) = 
\/£ 2 + to 2 we have A(£) < C Q A(£) 1_ I“I , we obtain 


C 


h(k,u)\<— Y l|A(Q) H - 1 V“ ff || il 


and 


C 


l“l<i 


\h(k,u)\<- Y HA(Q) |Q| - 1 V“ 5 || il . 


M<2 


Otherwise l 2 (k,u ) has support in a small cone around u, and we have 

I 2 {k,u) = f 
J R 3 


= f e^y^+^-^JuJ p k\u\)p 2 dp, 

J R+ 

Jv, P (k) — j 

Js 2 


with 


e lfc{ i - v -“ } g(pco)du;, 
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where g = X9- We obtain after an integration by parts 

r 2 7T 


Jv,p{k) — 7 


/ o 


,lfc{ 1 — COS' 


W} g(puj(e,(l))) 


dd 


0 


r /*7r 

+ kj 0 i e ' k{1 ~ C ° Smd *9(pu(0,4))dW0, 

Since we assumed g is supported in half cone around u, we have g(puj(9,7 r)) = 0. 
Hence we obtain 


Jv, P ( k )=r[ [ d c / > g(puj(9,(j)))d9 

to Jo Jo 

/»27T /*7T 

+ t / e lk ^ 1 ~ cos ^^d < j > g(pu>(6, (j>)) dc/)d0, 

to J o Jo 

and so 

MM) = irn-7 [ r r e Ife{v/ ^^” p|u|} ^(pa;(0>))#d0pdp 

\k\ M J R + Jo J-i T 

+ THTi f rr^-- ico - } ^(^))^^ ( 3 - § ) 

\k\ M /r+ do d-TT 

Let us now study the decay resulting from the dispersive behavior of the radial 
part. To this end, we follow the proof of the well-known Van Der Corput lemma. 
We study 

L(k,u, <t>, , 9) = [ e^Wp 2 +m 2 -PHco S ^)} d ^^ d ^ ))pdp _ 

J R+ 

Notice that, in view of (EH) , we are only interested by L(k,u,(f>,(f>,9) and 
L(k,u, 0,0, 9). First, for any differentiable function on R such that \ f\ > 1, 
we have for any a £ M + 


A({f £ R; |/(0I < a}) < a, 
for A the Lebesgue measure. We introduce 

h(p) = vV +m 2 - p\u\ cos ((/)'), 


(3.9) 


and we apply EUt to 


f(p) = 


dph(p) 


inf 

zesupp(g) 


{|0|L|k(l*l)l} 


We notice that dph(p) does not depend on u or (j/. With help of a smooth cut-off 
function, we split the integral L in two parts, one has support 

{p £ R+; \f{p)\ < a} 
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and the other is its complementary. In fact, we obtain exactly three interval 
corresponding to 

{p e R + ; f{p) < -a} , {p £ R + ; -a < f(p) < a} , {p £ R + ; a < f(p)} . 

In the first and third interval, we make an integration by parts and in the second 
interval, we use Estimate [EH . Hence we obtain the bound 

\L{k,u,<j),4>',0)\ < 


max < a, 


a\k\ inf 

KSsupp(g) 




max {p \d < j,g{pu}{6, </>))|} 
pSR+ 


We use 


and then choose 


a|fc| inf 

xesupp(g) 


f \d4,g{pui(9, 

0 ))l dp, [ p\d, 

R+ 

J R+ 

\ / 1 d^giruj 

(6», </>))| dr, [ r 

1./R+ 

J R+ 


1 


\f\k\ inf {|9| 2 x |/i(|x|)|} 

y xGsupp (g) l ) 


and plugging the resulting estimates for <// = 0 and <// = </> in we obtain 

estimate E3- 


3.2. Dispersive estimates for the free case equation. Thanks to the tools intro¬ 
duced in Section rm we are able to state the 


Theorem 3.1 (Dispersive estimates for free Dirac operator). For any p £ 

[1, 2], for all 9 £ [0; 1], for all s, s' £ R, such that s — s' > (| — 1)(2 + 9) and 
q £ [1, oo], we have 


— ltDrr 


lB i.^ BS , 


< (*(*))’ 


-l 


with 


and p' = -2 —■. 

1 p— 1 


K(t) = 


\t\- 1+9/2 if\t\ £ [0,1], 

i*r 1_8/2 ^ id £ [i,oo), 
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Proof. The proof is a straightforward adaptation of the one of Brenner in Ap¬ 
pendix 2 of IBre85| for the Klein-Gordon equation. We give a sketch of the proof 
for the reader’s convenience. Note that the proof of Theorem o for non zero 
potential, is based on the same ideas. 

We only need to prove the case p = 1, since the general case follows by interpo¬ 
lation of the case p = 1 and the charge conservation which corresponds to the 
case p = 2. Then using D m = V —A + m 2 ( 7T+ — 7r_) with n± = 1 K ±(D m ) = 
i {l ± \D m \~ 1 D m }, we obtain the estimates from those relative to the relativis¬ 
tic Schrodinger operator \J—A + to 2 : 




< m 


which in turn follow from 

Proposition 3.4. For any \ G ©(M^C 4 ), we define Xj( x ) = x(2 — |^c|)- Then 
for 9' £ [0,1], we have: 

1- if 0 ^ supp(x), 


|| e - lty /-A+m* Xj ^ v /_ Z \ + m 2^ \\ LljL „ < C- 2 ( 2 +«'b|i|-(l+e72) ( 3 . 10 ) 

where C is independent oft and j. 

2- if 0 € supp(x), 

|| e -nV-A+m* x{ _ A + m 2 )hi ^ Loo < C ( t )-3/2 (3.H) 


where C is independent t. 

We postpone the proof of Proposition IQ until the end of the proof of Theo¬ 
rem rrn 


We have 


e - lW -A+ m ‘ Xj(y /_ A+rn 2 )f 


<C^\\f\\l 


interpolating with Estimate of Proposition 13.41 for 9' = 0 when t < 1 and 

using Estimate lEHUl) with O’ = 9 for t > 1, one obtains 


2 JS 


—i ty/—A-\-r 


7 XjW~A + m 2 )f 


< C , K J -( t ) 2* , ||/|| 1 


with 


We use 


Kj(t.) 


2 j(2+e+ s '-s ) ( \t\ 1+8/2 if \t\ < 1, 
{it]- 1 - 0 / 2 if \t\ > 1. 


sup Kj < CK(t) 


if 2 + 9 < s — s' and Estimate iTTTl to prove Theorem ED Hence to conclude 
the proof, we need to give the 
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Proof (Proof of Proposition ^. f\j . Estimates of the same type, but for B \L P , L p j 

spaces with p £ [4/3, 2] can be found in lM^W80llMSW79l . In the present case 
p = la proof can be found in IB re 77j. This proof, which covers a much more 
general situation, is quite complicated. We propose here a simpler proof inspired 
by EMI. The kernel of + m 2 ) is given by ( 27 t) 3 ^ 2 Kj(x — 

y) where 

Kj(x,t) = f 
J R3 


Hence, we estimate the L°° norm of Kj. 


If |a:|/|i| < 2J- 1 /V2 2 T = T + m 2 or | x |/| 11 1, we use non stationary phase 

lemma in R 3 with help of the operator L = 


we obtain the estimate 


v/{2+m2 —ry • V. Hence, in this case, 


Vi' 2 + n 



e - 1 t v /?+mHx-f x .(^ 2 + m 2) d £ 


< C n 2~^ n ~ 3 ' > Pt\~ r ‘ 


for any n £ N. Otherwise, we apply Proposition 13.11 with h(r) = \Jr 2 + m 2 , 
k = t, u = x/t and g{x) = %j(|a:|). So if 0 ^ supp(x), to obtain 


f e-^^+^XjW? + m 2 ) dt, 
J R3 

lep '- 1 


c 

< — max 
" \t\ I/3|<1 


WxAVe + m 2 ) 


% i j 

R ± l€l 


df 


< 


C 2 2 -?' 

T 


Notice that in this case, \x\/\t\ > d > 0. If instead of Proposition 13.11 we use 
Proposition 13.31 with g = Xj, k = t and u = x/t , we prove the estimate 


\ K A*d)\ < 

The estimate ra is then obtained by interpolation. For imn . we use the 
classical stationary (Morse lemma) and non-stationary phase methods (integra¬ 
tion by parts) in R 3 . For more details about the method one can look at the end 
of the proof of Proposition 13.91 This ends the proof of 13.41 

This ends the proof of Theorem rm 


3.3. Distorted Plane Waves. Our aim is now to generalize the previous method 
to the perturbed case. Let us introduce the wave operators 


W ± 


lim e < Dm+v ) e - ltDm 


(3.12) 



































26 


Nabile Boussaid 


(for the existence and the completeness : Ran(VF ± ) = Ran(P c (£f)) of these 
operator see IGM01I Theorem 1.5]). With the intertwining property 

f{H)W ± =W ± f(D m ), (3.13) 

for any bounded borelian function /, and Fourier transform IF, we shall obtain 
for /i(£) = a • £ + to/3 

e ~' tH P C (H) = W ± e~ ltDm W ± * = W ± Te~ lth{Q) ( W ± T )* . 

So we can adapt the previous method if we are able to prove some estimates 
about the kernel i/>v of W ± T. The kernel ipv is called distorted plane wave. We 
notice that ipv is a 4 x 4 matrix valued function. 

We will show that the previous method works with ipv^vXj i n place of Xj with 
small modifications. So we need estimates on ipy Generally, distorted plane 
waves are studied like perturbations of free plane waves. So we will prove esti¬ 
mates on the perturbative part, written w in the sequel. 

3.3.1. Definition and properties. We need to introduce the free plane wave. Let 

h{k) = a ■ k + m/3 

for any k £ R 3 , notice that D m = h(P). This hermitian matrix has for eigenvec¬ 
tors the 

tpl(k,x) = e lk ' x u(k)ej 

where 

n(k) = (3.14) 

•\/2A(fc)(m + A (k)) 

with A (k) = Vk 2 + m 2 and (ej). are vectors of the canonical basis of C 4 . For 
more details see |Tha921 Section 1.4, Section l.F]. 

By definition, a distorted plane wave is a solution of the PDO equation 

(D m + V) ip = ± \/k 2 + m 2 if> (3.15) 

with for some j and any k G R 3 , i/j{k, x) — ^(k , x) tending to zero as x goes to 
infinity (in some sense), see m m75| section 5]. 

A solution of is a function ij){k, x ) of two variables here k is a 3-dimensional 

vector which is called the wave vector. A free plane wave i/j q satisfies the PDO 
equation in the case V = 0. Following |Agm75| , we introduce two families 

of function 

^v(k,x) = ip J 0 (k,x) - |i?+(A(/c))R(-)V , o( fc >’)} ( x ) 

for j e {1, 2} and 

v(k,x ) = ip J 0 (k,x) - {R^(-\(k))V(-)i/j J 0 (k,-)} (x) 

for j € {3,4}. The rest of the proof works also for Ry instead of Ry (the trace 
of the resolvent Ry was introduced in 112.111 1. 

In case there is no resonance at thresholds and no eigenvalue at thresholds, 
Theorem II. II gives us that Ry(X(p)) is in B(L 2 , L 2 _ a ) for any a > 5/2, this also 
work if a > 1 see Proposition 13.1(1 below. So the previous definition make sense 
if Assumption o holds and we have the 
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Proposition 3.5. Suppose that A ssum.ntinns i I . ?1 and i 1. & hold. Then for any k £ 
R 3 \ {0}, ijjy(k,x ) satisfies equation 13.1511 . 

Distorted plane waves define a generalized Fourier transform. We introduce 
ipv(k,x) £ Ad^C) the matrix with vector column ip y{k,x ) and we define 

C Fvf){k)=[ ipv{k,x)f(x) dx, 

Jr 3 


which is a priori defined on the Schwartz space <S(R 3 ,C 4 ) but will be extended 
to L 2 . 

Distorted plane waves are also called generalized eigenfunctions, since they cor¬ 
respond to “eigenvalues” associated with the continuous spectrum. Indeed, we 
can prove the 

Theorem 3.2 (Eigenfunction Expansion). The operator Ty defines a bounded 
linear map from L 2 into itself. Its kernel is given by the the sum of the eigenspaces 
of H. Moreover it is a unitary map from P C (H)L 2 onto L 2 (R 3 ) with 

(fyf)(x)= lim [ ip v (k, x)* f{k) dk, 

rwo ° Jk„ 

for any (An)neN a family of compact sets with K n C K n+ i and U ne NA' ra = R 3 . 
Finally, for any interval I Cl, one has 

Ui(H)f\\ 2 = [ \T v f\ 2 dk (3.16) 

J cr(/i(fc))n/7^0 

where a(h(k)) is the spectrum of h(k). 

Proof. The proof is an easy adaptation of the proof of |Agm75[ Theorem 6.2] 
(see also |R.S791 Theorem XI.41]), the main difference is that here we insert 
the unitary matrix u defined in iEHU) . Formula (I3T61) is nothing more than an 
adaptation of |Agm75| Formula (6.6)] or IR.S79I , Formula 82e’]. 

We also have the 

Lemma 3.2 (Intertwining Property). Let g be a bounded borelian function 
with support in R \ (—m, m) , we have 

T v g(H) = (goh)Tv (3.17) 

Proof. Using ism we obtain that (EH7I) is true for g = ij with I an interval 
of R \ (—to, m). We then obtain it for bounded borelian function with support 
in R \ (—to, m), usual density arguments and properties of functional calculus. 
More precisely, we use the fact that a bounded sequence of borelian functions 
which converges everywhere gives a sequence of bounded operators which con¬ 
verge strongly. 
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Hence we deduce that, for any \ £ Cg°(R), the kernel of e is given by 


[e ltH \(H)V C {H)\ (x,y) = [ ip v {k,x)*e lth{ - k) x(h(k))ip v {k,y) dk 



e~ ltH X (H)P c (H) = [T v )*e~ lth X {h)T v 


(3.18) 


We recall that we want to prove the decay of e ItH x(H) as t —> +oo in some 
Besov spaces. We observe that 


e~' th ^ X (h(k)) = e- ltx ^ X (X(k))P+(k)+e' tx ^x(-m)P-(k) 


where P+{k) (resp. P_(fc)) is the projector associated with the positive (resp. 
negative) part of the spectrum of h(k), i.e. 



Hence, in the following we study the functions 



(x, !/)el 3 xI 3 H / e TItA(fc) ( P±ipv{k , x))* (P±ip v (k, y)) x(K k )) dk. 


3.3.2. End of the proof of Theorem E3 We now prove Theorem ll.2l with help of 
three propositions which will be proven in Section [3 .3 .31 These propositions give 
some estimates on the perturbed part of the distorted plane wave. Following 
Cuccagna and Schirmer in ran . we write ipv{k,x) = e lk ' x {u{k) + uu(k,x)) 
where w is the perturbation part which satisfies 



(3.19) 


and we now state our propositions. 

Proposition 3.6. Suppose that Assumptions and fm hold. Then there ex¬ 


ists C > 0 such that for any k, x £ R 3 \ {0}, and any (3 £ N 3 with |/?| < 1, one 
has 



(3.20) 


Moreover one has 


|Vfctu(fc,a;)| < C 


(min{|a;|, |fc|}) 


(3.21) 


(k) (min{|x|, |fe|} |f|y — ]f[|) 
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We use this to prove the time decay in |t| 1 . Unfortunately this doesn’t work 
for the |f| -3 / 2 decay, hence we then study 


v(k, x)j 


e Ak\\x\ {ft+( +A (fc)) {Ve lk Qu(k)j}} {x), 
e —l|fe||ar| |i?+(_A(fc)) {Ve lk -Qu(k)j}} (x), 


if j £ {1, 2}, 
if j e {3, 4} 


and 


v(k , x)j 


e I l fe H x l |V fej R+(+A(fc)) {Ue lfc < 3u(/t)i}} (*), 

e _I ' fe " x ' {Vfci?y(—A(fc)) {Ue lfc < 3u(fc),}} (a), 


One has the 


if j £ {1, 2}, 
if i e {3, 4}. 

(3.23) 


Proposition 3.7. Suppose that A ssum.nti.on M 1. / 1 and \ 1 .‘A hold. Then if p > 3 + 
|/J| for some f3 £ N 3 , f/iere exists C > 0 such that for any k, x £ K 3 \ {0}, one 
has 


\V%v(k,x)\< 


C 


(1*1 


x k \ 

R \k\ ' 


Proposition 3.8. Suppose that Assumptions ^^ ] and \l.A hold. Then if p > 3 + 
|/3| /or some / £ N 3 , f/iere exists C > 0 such that for any k, x £ R 3 \ {0}, one 
has 


|v£«(fc,*)| < 


(min{|a:|, |fe|}) 


(fc) ^min{|a;|, |fc|} 


R 


R 


Using Proposition l3.(il 1X71 and ET%1 (which are proved in Section lX3.3l below). let 
us prove the following 

Proposition 3.9. Suppose that Assumptions \T~J\ and If.ill hold. Then we have 
for x £ Co°(K) with support in K. \ [—m; m] for any 9 £ [0,1] and j £ N, 


C9( 2 + e )i 

||e X(2 J H) || L i^ Loo < j^i+ 61/2 ’ (3.24) 

with C independent of t and j. 

We also have for x £ Cq°(M), f or an V 6 £ [0,1] 

\\e~ ltH x(H)P c (H)\\ Ll ^ < (3.25) 


with C independent oft. 

Proof. The proof works like the one of Proposition l3.4l with some modifications 
due to the fact that high derivatives in k of w(k,x ) grow with respect to x. 
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We need the L°° norm of the kernel of e ltH x(2 ^H). This kernel thanks 
to 13.1811 is given by 

Ij(t,x,y) = f {P + (0(u*®+w*(x, 0)} x 

Jr 3 

x {P + (0(u(0 + w(y, £)} x(2“ j A(0) dC 

+ / e +1 Vf 2 +™ 2 e -^-(*-!/) {P_(£)( U *(£) + «;*(*, f))} X 

Jr 3 

x {P_(0(«(0 +^,0}x(-2" J 'A(0)de- 

We notice that if we expand each integrand in terms of u and w, we obtain the 
sum of the integrals 

lf[z,z'](t,x,y) 

= [ e* lt Ve^e-'t<*-yHP±(0z%x,Z)P±(0z\y,0}x(±2- j m)dZ. 

Jr 3 

with z, z’ £ {u, w}. We notice that Ij~[u, x, y) + IJ [u, u](t, x , y) is the kernel 
of e~ ltDm Xj(D m ), hence we only treat the other integrals. 

For the Itl” 1 decay, if \x — y\/\t\ <C 2 J_1 / y/2 2 ^~ 2 + m 2 or \x — y\/\t\ 1, the 

phase has no critical point. We use an integration by parts in R 3 with help of 
the operator 



t , ^ 
y/e+m 2 


So with the estimate (EHU1) of Proposition 13. til and with 

dj * 

\Ji 2 + rrP 



we obtain the estimate 


\lf[z,z']{t,x,y)\<C2 2j \t\ \ 
with C independent of j and t. 

Otherwise if \x—y\/\t\ > c > 0, using first 13.411 of Proposition ^, ll and then 
of Proposition 13.(il we infer 

\lf[z,z'](t,x,y)\ < C'2 2j I tf” 1 , 

with C independent of j and t. 

For the |t| -3 / 2 decay, first if \x — y\/\t\ > c > 0, we write 

If[z,z'](t,x,y)= f , p \x_y\ )p 2 dp . 

J R+ 
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where 

J v (k) = [ e lk(1 -“' v) {P±(pu)z*(x, puj)P ± (paj)z'(y, pu)} x(±2~ J Hpu)) du. 

Js 2 

We can suppose v = (0; 0; 1) and so 


Jv(k) = iT £ eifc(i_c ° sw)) { p±(/9w(0 ’ puj{e ’ ^ ))x 

X -P±(pw(6»,0))z , (t/,pw(6»,(/>))|x(±2 _ ' 7 A(pw(6»,0))) sin (4>)d(j)d9. 
An integration by parts in <f> gives 


Jv{k) = J I"” [e lfe(1 ” COsW) { J P±(^(^^))^(^.P^(0,0))x 

x P±(pw(9, <t>))z’(y, pu{0, 0))}x(± 2_J A(pu;(0, 0)))] ” dd 
~ ^ jj £ g^d-co.(«)^{ P± ( pw(0 , .))*•(*, .))x 

x P±(paj(d,-))z'{y,puj(e, •))x(±2"' J A(pw(6», •)))}(</>) dcj)d9. 


The integrand of the hrst term can be rewritten in order to obtain a sum of 
two integral in <j> over the interval [0,7r]. To this end, we introduce a smooth 
cut-off function which splits [0,7r] in two parts one is a neighborhood of 0 and 
the other a neighborhood of 7r. Then most of the terms obtained after derivation 
can be treated by the method used for the |t| _1 decay. Only the two terms 
where derivatives of 0 and z' appear need a particular treatment. Now we have 
to distinguish the case z = z' = w from the two others where z = u or z' = u. 
If z = z' = w, the terms which need a particular treatment are bounded by Cltl^ 1 
times the supremum in ft of the e ([0,7r] x [0, 2 ir]) of 


</>,<!>')= [ e=F*tv^W-.p|x-,|co sW) {p^gp^x^u)}*. 
Jr+ 

X {P±{pw{9,<j>))d™z'(y, puj(9, </>))} x(±2 _J A {pu>(0,</>))) pdp, 
with n,m £ N such that n + m = 1. It is a sum of terms of the form 


f e it{ =F \J p 2 +m 2 -p ,x t yl (cos(0')-cos(0))-£iplfl +e i /pl|l} ^ 

Jr+ 

x \p±{pu)k,i X 

x {p ± (p W (M)k*' ; | x(±2- j A(p W (0,0))) pdp. 


where <j> is the angle between and the z—axis and ^(x, k) £ Ad 4 (C) is given 
by 

'(\x\\k\ + x ■ k)I& 0 C 2 

0 C 2 (-|x||fc| + x ■ k) J C 2 
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and Si, Si £ {±1}- We introduce 

K{p) = \p±{puj) kli (e- 1 'K k ' x tyz*(x,pu>)) ii } 

and 

| P ± (pu>(9, <f>)) ltk ' {d^z'(y,pu;(e,ct>))e^ k ^) ki J , 

Hp) = -FvV + ™ 2 ~ t ^ (cos (cj)) - cos (<j>)) - Sip y- + Ei'P-y 


and 


f(p) = 


d p (/)(p) 


inf {|9p0(]a:|,A)|} 

x£supp(xj) 


With help of a smooth cut-off function, we split the integral in two parts. One 
has support {t £ R; \f{t)\ < a} on which we use the estimate 

A({f £ R; \f(t)\ < a}) < a, 


for A the Lebesgue measure, since \f'\ > 1. The other is its complementary, in 
which we make an integration by parts. We obtain the estimate 


\jj~[r,r']{t,x,y)\ < Cot max {p\K (p)|| + , ] 2 x 

3 peA, at mf {\d 2 <t>{p)\} 

p£Aj 

x max | f {p \(d p K)(p)\} ; f \K{p)\ dp; 2~ J I {p\K(p)\} 

( •' A :i J Aj- j A j 

where A :j = g _1 {supp(xj)} with g(p) = y/p 2 + to 2 . Hence with (1H.21D of Propo¬ 
sition ITiil Proposition EJ and decay of derivatives of P± , we can choose a = 
2 2 i \ft and we obtain the bound of in this case. 

For the case ( z,z') = ( u,w ) (the case ( z,z') = ( w,u ) is similar), we study by 
the same way the integral 


I e lt { : Fv / P 2 + m2 -P^T^( cos ( , ? i, )-cos(^))-e i /pl|l} x 
/R+ 

X | {p±(pw)fc,i (d^z*{x,pu})). l j x 

x jp ± (pu;(^))^ (d^z'{y,pu:(e^))e^ k ^) ki X x 


x X(±2 J X(pu(6, (/))))} pdp, 


where </> is the angle between ^ and the z— axis. 

If \x — y\/\t\ <C 1, we can suppose |x — y\/\t\ < |£|/(2-\/£ 2 + m 2 ) for any £ £ 
supp(xj) and instead of applying the trick of the proof of Lemma mi f integration 
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by parts with respect to an angular variables) to the integral Ij[z,z'](t,x,y), 
we make an integration by parts with help of 

%l 

isl _|_ J_ . x-y ' 
y^+rn 2 l?l 4 

The rest of the proof is the same. 

We now turn to the proof of estimate i;-).25ii . the kernel of the operator is given 
by a sum of terms of the form 

lf[z,z'](t,x,y) = f (*-!'>x 

J R 3 

x{p ± (Oz*(x,op±(ty(y,<;)}x(±m)dt. 

We first notice that Proposition Id.(il implies that this integral is bounded. Then 
we split the integral in two parts. One is supported in a small neighborhood of 
the critical point of the phase, the other is its complementary. To treat this last 
integral we work exactly like the case “ ' x ~ v ' <C 1”, just mentioned above. For 
the other one, we apply the Morse lemma to reduced the study to 

[ e^« 2 {p ± (f(o)z*(x,m)p±m))Av,fm x(/(0 )% = 

J R 3 

/ / pe Tltp2 [p±(f(puj))z*{xJ(puj))x 

Js 2 J R+ L 

X p± {f(pu))z'(y, /(pw))|x(/(pw)) dpdu, 

where \is the product of an indicator of a small neighborhood of the critical 
point with x(±A(-)). Then an integration by parts in p and the Van Der Corput 
lemma give (P5l) when 9=1. Since we have that the integral is bounded the 
general case easily follows. 

We are now able to write the proof of Theorem PI using Proposition |^J 
Proof (Proof of Th.eorem \l.2\) . We notice that 

= D-^kiDmW^H) - D-'cfkiD^VfaiH) 

= 2~ k MD m )2 j fa{H) - 2~ k 4 )k {D m )V(f 0 (H) 

We can also use since the support of (f> 3 is far from 0 

MDmWjiH) = - </> k (H)VH- 1 ^(H) 

or higher power in -D” 1 or H _1 to obtain with (Id. 2 511 

/^o( 2 +6)j 

< 02-^-* tl+g 2~ r ' J - fc ' 

for any reals r, r' with C independent of i, j. Hence if r, r' > 0, we work like 
in the proof of Theorem rOlfic like in Appendix 2 of }Hre77j l to conclude the 
proof. 

It now remains to prove Proposition Id. 61 Id. 71 and Id. 81 
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3.3.3. Some estimates. 

Estimates for w. We remind us of the definition of w in cnsu and we introduce 
R$(k) = e~ ikQ R+(±X(k))e ikQ . (3.26) 

We have 

Lemma 3.3. For any a € N 3 , let be o > 4 + |a|. Then there exists C > 0 such 
that for any k, x £ R 3 \ { 0} 


|{v£i? o ± (fc)(Q)-g}0r) 


< 


C 




Ja| 


x _ k 

FI FI 


|W 2 +I“I> 


< fc > H < | X | 

We also have that there exists C > 0 such that for any k, x £ I 3 \ { 0} 

\{v%R±(k)(Q)-°q} (x) 


(3.27) 


< C 


[x) 


Ot —1 


(min{|a:|, \k\}) 


( fc) a_1 


(k) (min{|x|, |fc|} 


FT FI 


lvv 2 +i“i, °° • (3.28) 


Proof. We write 


(f? 0 ± (fc)(O)- g ) (x) 

f e 1 ^ ± Fllyl+ fe '^} ( a ■ (x — y)q{x — y) a-'Vq(x—y) 


4tt|2/| 


(x - y) a + 2 


+ (a ■ k + to/ 3 ± A(fc)) f 

Jr 3 


(x - y)° 

a i{±\k\\v\+k-y} 


dy 


vi q(x - y) 


A-k\v\ ( x-y) a 


dy. 


We restrict our study to 3?(j"(fc) since the two cases are similar. Hence we only 
need to estimate integrals of the form 


R(k)(x) = f e 1 

J R3 


{\k\\y\+k-y} 


- y) 

\y\ 


dy 


with u G W’i +I “ I ’ 00 (M 3 , C). 

In a first step, a straightforward calculation shows that 
|V^(fc)(x)|<C(x)H~i|| u || 


(3.29) 


if a > 3 + max{|a| — 1; 0}. Then using the trick we used in the proof of Propo¬ 
sition EU we obtain 


V a k R(k)(x) = 


JfcjW 


pdFIFI +k-v} 


1 R 3 


(viQir 


u(x — •) 

~w~ 


(y) dy , 
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and so with (l3~29l) . we infer 

|V£i2(A0(z)| < (3.30) 

since a > 3 + max{|a| — 1, 0}. 

In a second step, we apply Estimate (1331) of Proposition 1331 to R(k)(x), this 
gives 


|Vgi2(fc)(*)| 

C 

< 


max 


|fc|M + 1 |/3|<1 + M ] J R 3 


I y\ 


1/31-1. 


k _ y_ 

R \v\ 




Hence we need to estimate on integral of the form 

1 1 1 


G(x,uj) = [ 

Jr 


R 3 \y\ n (x - y) s 


Jy\ 


— to 


- ■) 

\Q\ 


dy 


(y) 


dy } . 


with ui £ S , — |or| + 1 < n < 2 and s > <r. To obtain appropriate estimates, we 
use 


1 


_x _ y_ 

M \y\ 

to write for 9, 9' > 0 such that 9 + 9' = 1, 


\x~y\> - max{|y|, |x|} 


1 


G( *’“ )GC Lww H 


\0s 


x _ y_ 

R \v\ 


dy 


< 


< 



C 

f 1 1 

M (m 

UJ — y -T 

^ Jr+ r n 2 (|a/| ~r) es 

C(x ) |a|+1 


\ x \ (M 

(jJ — jy 

M 

y 


dr 


if 9's > 2 and 9s > 1 + max{2 — n: 0}. Since G(0, u>) is bounded, we obtain 

C{x) |q| 


G(x,<jj) < 


Hence, we obtain with estimate (TOol) 


C 


(\ x \ 

, , _ X 

\ 

V 1 

M 

/ 




\i«i 


|fc| \x\ 




(3.31) 


which gives estimate (13371) . 

In a third step, if k/\k\ ^ x/\x\, we split the integral for \7%R(k)(x) in two 
parts with help of a smooth cut-off function defined in S 2 the support of which 
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is a half cone determined by the bisector plane of the couple (k/\k\\x/\x\). So 
we obtain V%R(k)(x) = R\(k)(x) + R 2 (k)(x) with Ri(k)(x) having a support 
containing x/\x\ and R 2 (k)(x) having a support containing k/\k\. We then apply 
the estimate (TTfil) of Proposition EP to Ri(k)(x) to obtain 


I V%Ri(k)(x)\ < 


c 


H+2 


k x 

|fe| R 


x max 


\0\<2+\a\ | J R 3 
This gives the estimate 


I y\ 


1 / 31 - 2 . 


k _ y_ 

l fc l \v\ 


V' 3 


u(x — •) 


(y) 


dy ) . 


C (x)l Q l -1 

|V fc i?i(fc)(a;)| < |^|| q | +2 l 2 ll u llwi“ l+2 ’' 


|fc| \x\ 


since a > 2 + |ct|. Using iron . we infer 


c 


( fc ) H + 1 (y/\k\\x\ 

k _ x 

R R 

) 2 




or, using iliHUli . 


I V%Ri(k)(x)\ < 


C 




-1 


( fc ) H <|*| 


k x 

R R 


\\u\\ w m+2,™- 


For R 2 (x)(k), we use the inequality \x — y\ > 1 ^ Tl to obtain 

C(x)l a l- 1 


R 2 {k)(x) < 


(fc)l Q l(|a 


|fe| |aj| 


\ u \\wl^ 


since a > 3 + s + max{|a| — 1, 0}. So now we easily deduce estimate (13.281) . 
For the sequel, we need the following 

Lemma 3.4. Let be s £ R and <fi a C°° function such that there is a > 0 with 

\/a G N 3 , |V“0(*)| < 


We have that [( P ) s , <f{Q)\ is bounded from H * into H*, with q'+a > q and t'+l > 
t + s. 
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Proof. We want to prove that 

(Q) q (P) t [(P) s ,HQ)}(P)- t '(Q)- q ' (3.32) 

is bounded in B(L 2 ). Using the identity 

[(py^m = i(p) s / 2 , <t>(Q)\(p) s/2 + (p) s/2 [(p) s/2 ,<t>(Q)] 

we reduce the proof to the case |s| < 1. And with the identity 

i(py,cfm = -(pn(p)- s Am(p) s 


we only need to study the case — 1 < s < 0. The proof in this case is based on 
the following identity for — 1 < s < 0 


<p>- = (-/i + ir' 2 = ^d!Lill) 


+°° { 

vn 2 j 


—A + 1 + w 


dw. 


So we have 


[<p) s , <t>(Q)\ = jr +l ) s/2 - kAd -A + i{<t>m +R-n 


with 


Rm — 


(—l) m sin(7r {§}) [ + °° 


{-A + l + w) m+1 


Ady+lMQ)) 


dw 


— A + 1 + w 


Then we use = 1 — , and we commute powers of (P) with 

operators of the form S7 a <p(Q). Hence we can repeat the previous computation 
until we obtain only non positive powers of (P) in (Id.d2ll . So we only need to 
prove that operators of the form 


[{Q)\<f{P)]{Q)~ q ' 

with q < q' + 1 and </> satisfying the assumption of the lemma are bounded 
in B{L 2 ), we just repeat the previous calculation but we switch the role of P 
and Q. This ends the proof. 


We now state a particular version of the Limiting Absorption Principle for H. 

Proposition 3.10. We assume that Assumptions I7.il and ro hold. Then for 
any a > 1 there exists C > 0 such that for any k £ K 3 

\\RvWh(Ll,L>_„) < c. 

Proof. In fact, we just need to prove that for any a > 1 there exists C > 0 such 
that for any A € R \ (—m, m) 

ll-Ry(^)llB(LUPLJ — C- 

Using Theorem 11 .1 1 we have that it is true if o > 5/2. Then we use Born 
expansion 

P+(A) = P+(A) - P+(A)VP 0 +(A) + P+(A)VP+(A)VP+(A) 
and IIM99 s Theorem 2.1 (i)] to end the proof. 
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We are now able to give the 

Proof (Proof of Proposition \3. 61) . We only give a the general idea of the proof 
and we leave the details to the reader. We notice that with Ry defined by (13.261) . 
we obtain 

w = RyVu 

with an abuse of notation since we avoid to distinguish the case where we 
have Ry or Ry. We recall the identities 

RyV = R 0 V - RqVRvV = R 0 V - RqVRq + R 0 VR v VR 0 V. (3.33) 

Since, we have 

Rv = (1 + RoVy'Ro, (1 + RoV)- 1 = 1 - R v V, 

for | a | = 1, we obtain 

V c (R v = (1 - R v V)VlR 0 {l - VR v ). (3.34) 

Using l ETSfl . we obtain a formula where only derivatives of Rq appear (if there 
is derivatives). Then between a derivative of Rq and Ry, we insert a Rq with 
the identity lliOHll : 

RvVV^RqV = RoVV^RqV - RqVRqVS/^RqV + R 0 VR v VRoVX7%R 0 V. 

This ensures that if p > 5, V or its derivatives decays enough to use Esti¬ 
mate iETTHi and Proposition mm Since these estimates need derivatives and 
Sobolev’s injections, we apply Lemma, 13.41 to conclude the proof. 


Estimates for v. We remind us of the definition of v in El and we introduce 

S(f’ E2 {k) = e- Eie2l|fc|IQI i?^(e 2 A(fc))e lfc ' Q , 

where gy € {—1, 1}- With an abuse of notation, we will write v = SyVu. We 
have 

Lemma 3.5. There exists C > 0, such that for any k £ R 3 \ {0} and (3 € N 3 

C 


(vf5 0 £l ’ £2 (fc)(Q)-^) (x) 


< 


x k 

M 16 


I w 2 + 101. 


for any a > 3 + \(3\. 

Proof. 

(S £l ’ £2 (fc)(Q)- g ) (*) = 

el e 1 e 2 {\k\\x-y\-\k\\x\+e 1 e 2 k.y} f a . yq ( y j a . > 


/ R 3 


47r|a; — y\ 

(a ■ k + m/3 ± A (k)) ( 

J R 


(■ y ) a+2 ' ( y)° J dy 

e iEie2{|fc||*-y|-|fc||*|+eiE2fc-y} q l y ^ 


4:n\x — y | 


(■vY 


dy. 
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For the sake of simplicity, we only write the proof when (3 = 0. The proof for 
derivatives works in the same way using ||x — y\ — |x|| < \y\ and a > 3 + \(3\. 
But the proof for the case (3 = 0, has been already done since 

| [Sl^{k){Q)-°q) (x)| = | (R s 0 ^(k)(Q)-°q) (x)| . 

Hence using Proposition 13.101 we able to write the 

Proof (Proof of Proposition |3771 ) . We write with an abuse of notation 

v = SyVu, 


and we use the Born formula 

S v V = S 0 V - S 0 VR v V, 

together with Lemma Id. 51 Propositions EH1 and Ol The proof works like the 
one for w. 


Estimates for v. We remind us of the definition of v in (I3T771) and we introduce 

Ty 1 ’ 62 (fc) = e- Eie2l ^ Ql+lkQ V k R^{e 2 X{k)) elk - Q , 

where € { — 1, 1}- With another abuse of notation, here we will write v = 
TyVu. We have 

Lemma 3.6. There exists C > 0, such that for any k £ R 3 \ {0} and (3 £ N 3 

(min{|x|, \k\}) 


(v 0 k Tt £2 (k)(Q)-°q) (:X ) 


< C- 


(k) ^min{|x|, |fc|} 


x k 

l*l i fe i 


2 II9II W 2 +I^l, °° s 


for any a > 4 + 


Proof. This is an obvious adaptation of the proof of Lemma 13.51 we just notice 
that one has 

| (Tt e2 (k)(Q)~ a q) (z)| = | (Vfc-Ry (e 2 A(fc))(Q) _CT g^ (x) 


Hence, we have 

Proof (Proof of Proposition Id.#)) . One more time, we write with an abuse of 
notation 

v = TyVu + S v VX7 k u, 

The second term of the right hand side could be studied exactly as we done in 
proof of Proposition 13 . 71 and for the first one we use the formula 

T v V = T 0 V - ToVRyV + S 0 W k RyV, 

together with Lemma lib 61 Propositions m and rrrui The proof works like the 
one for w. 
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4. The linearized operator 

In this section, we study the spectral properties of the linearized operator, asso¬ 
ciated with Equation o> around a stationary state. This will be useful since 
we compare the dynamics associated with Equation Q to the dynamic of the 
linear Dirac equation associated with H. This comparison is possible only be¬ 
cause when the PLS is small, the linearized operator is a small perturbation 
of H. 


4-1- The manifold of the particle like solutions. First we notice that Proposi¬ 
tion o which gives the existence of stationary states, is a consequence of 

Proposition 4.1. Let H be a self adjoint operator on P 2 (R 3 ,C 4 ) and with a 
simple eigenvalue Ao associated with a normalized eigenvector <j)Q. Assume that 
there is a neighborhood O C R of Xq such that for all A G O the operator (H — 
A)- 1 P 0 is mP(P 2 (R 3 , C 4 )) for any o G R+, and in P(P Z (R 3 , C 4 ), P /+ 1 (R 3 , C 4 )) 
for any l G N, where Pq is the projector into the orthogonal space of <j> o. Let F G 
C fc+ 1 (C 4 ,C 4 ) such that F(z) = 0(\z\ 3 ). 

Then for any a G R + , there exists fl a neighborhood of 0 G C, a C k map 

h-.n^ {<M X n h 2 (R 3 , c 4 ) n (R 3 , c 4 ) 

and a C k map E : fl i— > R such that S(u) = u(f> o + h(u) satisfy for all u G 12, 

HS{u) + VF(S(u)) = E{u)S{u ), 

with the following properties 

f h{e i6 u) = e iB h(u), V6> G R, 

I h(u)=0(\u\ 2 ), 

| E{u) = E(\u \), 
l E(u) = A 0 + 0(|u| 2 ). 

The proof of this proposition is an obvious adaptation of the one of iPW97l 
Proposition 2.2], and we don’t repeat it here. One can also obtain it by means 
of the Crandall-Rabinowitz theorem but it doesn’t give immediately the decom¬ 
position associated to the spectrum of H = D m + V. 

To show that {H — A) _1 Po is in £>(P 2 (R 3 , C 4 )) for any a > 0 , we just need to 
prove that a i—>■ e“^ (H — A) - 1 Poe -a ^ is of class C k near 0 in B(L 2 ( R 3 , C 4 )) 
for any k G N, this can be proved with help of of [HisOOl Lemma 5.1]. To prove 
that (H - A) _ 1 P 0 for any l G N is in B(H l ( R 3 , C 4 ), P* + 1 (R 3 , C 4 )) for any l G N, 
we first notice that (D m — A ) -1 is in B(H l ( R 3 ,C 4 ),P Z+ 1 (R 3 ,C 4 )) then we use 
wave operator, see Id. 121 and IGM01I Theorem 1.5], and the intertwining property, 
see l3~T3l to conclude. 

We shall need some properties of stationary solutions of (11.31) . Following lEtisOOl , 
we have the 

Lemma 4.1 (exponential decay). For all f3 G N 2 , s G R + and p,q G [1, oo]. 
There is 7 > 0, e > 0 and C > 0 such that for all u G Bc(0,e) one has 

\\e^d^S(u)\\ B s q <C\\S(u)\\ 2 , 

where = d u l3iud u 2 Q u ■ 
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Proof. In fact we prove that for any k in N there is 7 > 0 and £ > 0 and C > 0 
such that for all u € Bc(0,e) one has 

||ert«>a£S(ii)|| ff * < C||S(u)||2. 

Then interpolation and the following property of Besov spaces over R 3 permit 
to conclude: B 22 = H s , B a r C B qi ' q if s' < s, B‘‘: q C B a q ifl<r<p<oo 
and u-n/r = s - n/p and \\uv\\ B ^ q < if | + § > \ + y- 

We only prove the lemma for (3 = 0, the other cases are similar. We have 

D m S(u) + VS(u) + VF(S(u)) = E(u)S(u). 


Let us introduce the M—linear operator W of multiplication by the matrix val¬ 
ued function x £ R 3 1 —> —iDVF(S(u)(x))i + V(x). We obtain, with the gauge 
invariance of F, the identity 

WS(u) = VF(5(u)) + VS(u). 


The “potential” W tends to zero as x goes to 00 . In fact, as a function of x, W 
is in L 1 n L°°; we can write W = W c + Ws where W c is compactly supported 
and \\W s \\ L i nL oo < 6. 

We have that D m + W ,5 — E(u) is invertible for 5 sufficiently small and 
e^ Q) S(u) = e 7<Q> (D m + W S - £(u )) _1 e - 7 < Q) {e 7(Q> W c S(u)}. 

For 7 small, ^ D m + + Ws — E(u)j is invertible in L 2 and 

+ W s - E(u)j e 1 ^' l W c S{u). 

This proves the lemma for k = 0 since e 7 '® W c is bounded. Now we notice that 


e l{Q) S(u) = + 


|P| ( Dm+I l^ + WS ~ E{U) ) 

= £ - £ £ w + w ‘ - E( “») ( D " + ’¥ + m ~ EW ) 

Hence we obtain 

|| e 7 <Q>5( u )|| ii . fc < C'||5(u)|| ff k-i. 

This identity proves the lemma by induction. 
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f.2. The spectrum of the linearized operator. Here we study the spectrum of 
the linearized operator associated with Equation JE3 around a stationary state 
S(u). Let us introduce 

H(u) = H + d 2 F{S(u )) - E(u) 

where d 2 F is the differential of VF. The operator H{u) is R—linear but not 
C—linear. Replacing F 2 (R 3 ,C 4 ) by L 2 (R 3 ,R 4 x R 4 ) with the inner product 
obtained by taking the real part of the inner product of L 2 (R 3 ,C 4 ), we ob¬ 
tain a symmetric operator. We then complexify this real Hilbert space and ob¬ 
tain F 2 (R 3 ,C 4 x C 4 ) with its canonical hermitian product. This process trans¬ 
forms the operator — 1 into 

7 _ ( 0 -Id C i\ 

J -\ld C 4 0 J ■ 

For 4> £ F 2 (R 3 ,R 4 x R 4 ) C F 2 (R 3 , C 4 x C 4 ), we still write <f> instead of 

\^J- 

The extension of H{u) over F 2 (R 3 ,C 4 x C 4 ) is also written H(u) and is now a 
real operator. 

The linearized operator associated with Equation Q around the stationary 
state S(u) is given by JH{u). We shall now study its spectrum. Differentiat¬ 
ing Id, we have that 

M “ =Spa ”{al; s< " ) ’ ds s(o) } 

is invariant under the action of JH(u). We notice (see l(!NT04l 'l that 

Ho(u) = Span { JS(u), 9| u |5(u)} . 

Using gauge invariance and differentiating, we obtain 

JH(u)JS(u) = 0 and JH(u)d\ u \S(u) = d\ u \E(u) JS(u). 

Hence Hq(u) is contained in the geometric null space of JH{u ), in fact it is 
exactly the geometric null space as proved in the sequel of this subsection. First, 
we see that JH(u) has two other simple eigenvalues, as stated in the following 

Lemma 4.2. Let be 


st( 0) 



and S 1 (0) 



Suppose that Assumptions \ 1 . 1 \\ 1 . 2 \ hold, then there are e > 0 and four C°° 
maps Ef : £?c(0, e) t—> C and kf : Bc(0,e) t—> {S'j t (0)} ± such that 

JH(u)Sf(u) = Ef(u)S±(u), 


with || S^~ (w) || = 1, 

St{u) = Sf(0) + kf{u), 
with Ef(u) = ±*(Ai — Aq) + 0(|it| 2 ) and kf( 0) = 0. 
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Proof. This can be proved in the same fashion as [PW971 Proposition 2.2] using 
Assumption 11,31 

We also obtain 

Lemma 4.3 (exponential decay in Besov spaces). Suppose that Assump- 
tions I7~71| i ,J\ hold, then for any (3 £ N 2 , s £ R and p,q £ [l,oo] there 
is 7 > 0, £ > 0 and a positive constant C such that for all u £ Bc{0, e), 

\\eiMdZst(u)\\ B s q < q|S±(«)|| 2 , 


where d u = d^Jud^u - 

Proof. The proof is exactly the same as the one of Lemma rm 

Let TL±\{u) be the space spanned by Sf (u). Let us now prove that the orthogonal 
space with respect to the lrermitian product associated to J 

H c {u) = {Hq(u) © Ti+i (u) ® 

contains no eigenvector. We notice that R c (u) is invariant under the action 
of JH(u). We have 

Lemma 4.4 (Continuous subspace property). If A ssumvtions \T7\ \l.j\ hold, 

let P c (m) be the orthogonal projector onto Tt c (u). Then there exists e > 0 such 
that for u', u £ £?c( 0 , e) 


p c({u))\ Hc{ul) : H c (u') i * TLciu) 

is an isomorphism from B^ q { IR 3 ,C 8 ) D TL c {u’) into B^^M 3 , C 8 ) ft H c (u), for 
any s £ K. + and any p,q £ [l,oo]. The inverse R{u',u) is continuous with 
respect to u and u'. 

Proof. This proof is a straightforward adaptation of the one of |CNT041 Lemma 
2 . 2 ].' 


So we have 

Lemma 4.5. Under the assumptions of Provosition \T~T\ there exists e > 0 such 
that for any u £ B c( 0 , e) 

\\{Q)-°e° JH ^P c {u)n < C(s )- min ^ 3 / 2 }||<Q)>||, € Ll 

[ ||(g)- <r e 8jff (“>P c (u)^ir ds < C\m, Vip £ L 2 . 

JR 

As a consequence, TL c {u) does not contain any eigenvector. 

Proof. For the sake of clarity, we introduce 

c(u) = J JL s{u) ’ J5i+(m) ’ jsr(u) ) ■ 
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Writing Duhamel’s formula for H(u) with respect to H — E(u ), we obtain 
e tjjf(«0p c ( u ) = e w(fl-£W)p c ( u ) 

+ [ e {t - s)J(Dm+v ~ E{u)) Jd 2 F(S{u))e sJH{u) V c {u)ds. 

Jo 


We have 


P C (0)Ih c („o = £(«', 0) = Id L 2 + £ Mu', 0)) <0(0)1 

i 

where the coordinates of ai(u',u) are linear combination of the coordinates 
of <(u), so it can be extended to L 2 _ a and we have 

\\(Q)-"e- tJH Mp c (uW\\ 

< IIW Pc(0)|^ ( „,) (Q)"||{||(O)- ff Pc(0) e - tJ ( D - +y - B (“»P c ( M )^|| 

+ 

< c , (t)- min{<T ’ 3/2} || (<3)^11 

+C f (t — s )- min l CT ’ \\{Q) 2cr D~VF(S(u))\\\\(Q)~ ,J e~ isH( ' u ^P c (u)iIj\\ ds 

Jo 

We then introduce 

M(t) = sup {(s)- min{CT ’ 3/2} ||(0)- <T e- sJff(M) P c (w)V’||} 
se[o,t] 

and we obtain for \z\ < e 



M(t)<C(\\{Q)"i/>\\+eM(t)) 

which gives for e sufficiently small 

M(t)<C\\(Q)^\\, 

or 

\\(Q)~ a e~ sJH( - u) V c {u)if\\ < C(s}~ min l cr ’ 3 / 2 l||<Q) < 7 V>||. 

With the same method, see Lemma we obtain the second estimate. 

Then we obtain with the second estimate that there is no stationary state in 
the range of P c (it) that is to say H c (u). 

This gives 

Lemma 4.6. We have, for sufficiently small u £ C, E e (u) € iR with E±(u) = 
—E^(u) and Sf (u) = S E (u) for the conjugation of C 8 . 
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Proof. The last statement straightforwardly follows from 

JH(u)Sf (u) = E±(u)Sf(u), 

since there is no more eigenvalues than the 0 and E±(u), we obtain E±(u) = 
Ef(u). 

Then we specify the essential spectrum of JH{u). A classical study gives that 
the continuous spectrum of JH( 0) is given by 

{iA; A £ R, |A| > min{|m — Ao|, | m + Ao|}} . 

Using Weyl’s criterion (see |RS78I Theorem XIII. 14, Corollary 1], the adaptation 
is quite easy in our case), we obtain that the essential spectrum is 

{iA; A £ R, |A| > min{|m — E\, \m + E\}} . 

Hence Ef(u) are necessarily purely imaginary. Indeed if H(u) — Ef(u)J is 
not invertible then H(u ) + FTf (u)J is not invertible too. Since — E±(u) is not 
in the essential spectrum, it is necessarily an eigenvalue in the neighborhood 
of ±i (Ai — Aq). Hence this gives — Ef-(u) = Ef(u). 


4-3. Decomposition of the system. We want to decompose a solution <f of the 
equation Oil with respect to the spectrum of JH(u). And in fact, we only 
study the resulting equations for these different parts of the decomposition. 
First we isolate a part which corresponds to a PLS. For any solution of 
over an interval of time I containing 0, we write for t £ I 

(fit) = e- 1 ^ ds (, S(u(t )) + v {t)). 

In order to give an equation for rj, we introduce the following space 

K(u) = j»j £ L 2 (R 3 , C 8 ), ^ Jr,, = 0, (jp, -^S{u)j 

In fact it is the space 

H+i(u) © 7t~i(u) © H c (u) 

which is invariant under the action of JH{u) and we state the 

Lemma 4.7 (decomposition lemma). Let be s > 0 and p > 1 there exist S > 
0 and a C°° map U : Bw s ’p{ 0, S) e-> Hc(0, e) which satisfies for if £ Bw s -p( 0, 6) 

if = S(u) + r], with r\ £ Hq(u) u = U(tf). 



Proof. It is ICNT04I Lemma 2.3]. 
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This lemma ensures that we can impose the orthogonality condition 

V(t) G (4-1) 

So instead of solving the Equation d in (/>, we want to solve the equation 

id t r] = {H - E(u)} rj + {VF(S(u) + 77 ) - VF(S(«))} - idS(u)u . 

= {H+ d 2 F(S(u)) - E(u)} r) + N(u, 77) ~ i dS(u)ii [ > 

for 77 £ 7ij~(u(t)). Here d 2 F is the differential of VF and dS the differential of S 
in JR 2 . To close the system, we need an equation for u. Let us now derive an 
equation for the path u, by means of G31>: 

(' ri(t),JdS(u{t ))) = 0. 

After a time derivation, we obtain 

0 = (J H (u(t))r](t) + JN(u(t)),r](t )) 

+ dS(u(t))u(t), JdS(u(t))) — (77, Jd 2 S(u(t))u(t)). 
Since S(u) £ JTio(u), we have 

(H (77)77, dS(u)) = (77, H(u)dS(u)) = (77, dE(u)S(u)) =0, 

we obtain 

[(JdS(u(t)),dS(u(t))) - (. Jrj(t),d 2 S(u(t)))]u(t ) = -(N(u(t),rf{t)) : dS(u(t))). 
So we notice that 

[(JdS(u(t)),dS(u(t))) - (Jrj(t),d 2 S(u(t)))] = (^ ^ + 0(\u(t)\ + ||r?(i)|| 2 ), 

which proves that [( JdS(u(t)), dS(u(t))) — (Jrj(t), d 2 S(u(t)))} is invertible for 
small \u(t)\ and ||^(t)|| 2 , we therefore introduce its inverse 

A(u,rj) = [( JdS(u), dS(u)) — (Jr], d 2 S'(Tt ))] -1 

and write 

d t u(t) = -A(u(t),r](t))(N(u(t),r](t)),dS(u(t))). 

Plugging in Equation and similarly to the linear case we decompose 77 with 
respect to the spectral decomposition of H(u) = H + DX7F(S(u)) — E(u) 

r](t) = a + (t)Si (u) + a~(t)Sj ( u ) + z(t ) 

with 2 £ Ti c (u) n F 2 (R 3 ,R 8 ) and a~ = a + . We obtain the system 

'u =—A(u,i])(N(u,r]),dS(u)) 
aA = -E ± (u)q ± + (JN(u,rj), JS±(uj) 

+(dS(u)A(u, r))(N(u, 77), dS(u)) JS±(u)) 

< -((dSf(u))A(u,r])(N(u,ri),dS(u)), JS'i t (u))a ± 

-{(dSf (u))A(u,ri)(N(u,ri),dS(u)), JS^(u))a T 
d t z = JH(u)z + P c (u) JN(u, 77 ) 

+P c (u)dS (u)A(u, r]){N(u, rj), dS(uj) 

~(DP c (u))A(u, ij)(N(u, rj), dS(u))rj 

which we will now study. We notice that this equation is defined only for 2 small 
with real values, a~ = a + small and u small. 
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5. The stabilization towards the PLS manifold 


We now build a solution which stabilizes towards the manifold of the stationary 
states. To this end, we will use Theorem 11.11 and Theorem o to prove that z 
tends to zero in L°° and Lf oc . It is possible here since we build solutions for 
which we ensure that a + and or also tend to zero. We do not think that this 
convergence holds for all initial states but we do not know any counterexample. 
We also notice that we look for a real solution <j> = S(u) + ij, hence rj should be 
real and therefore a~ = a + . 

We impose the following condition 


C 

' '-W' 


Under the assumptions of Theorem Ol let us define for any e, 6 > 0 


U(e,5) = lu £ C°°(R, Bc(0, e)), lim u(t) = r«oo exists, 

t—► + OO 


S 2 

I U{t) - Moo| < ^,Vf > 0 


and for any u £ U(e), let s, s', (3 be such that s' >s + 3>/3 + 6 and a > 5/2, 
we define 


Z(u,5) = {z£ C°°(R, L 


OO/TU) r2/'TIj3 IP8 


)), z(t) £ H c (u(t)), 


max 


sup sup {{v) 3 / 2 \\z(v)\\ b p }, 


i>G[0,+oo] 


vG [0,+oo] 


sup {(w) 3 / 2 {||z(V)||ir 1 J 

fG[0,+oo] 


<sy 


Then we define the space 


f2(6) = fa = (a + , a ) G C°°(M), a = a+, sup (f) 3 / 2 |a(t)| < 5 2 1 . 
I teR+ J 


5.1. Step 1: Construction of a. For any u £U{e, 5) and z G Z(u, S), let us define 
a map G u ,z on f2(5) by 

GuA^it) - - 

+ (dS(u(v))A(u(v),r](v))(N(u(v),ri(v)),dS(u(v))) St{u{v))) 

- ((dSt{u{v)))A(u(v),r]{v))(N(u(v),r](v)),dS{u(v))), Si(u( v ))) a± ( v ) 

- ({dSf(u(v)))A(u(v),r](v))(N(u(v),r](v)),dS(u(v))), S^(u(v)))a T (v) 1 dv. 
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We want to show that Q UyZ stabilizes 12(5) and is a contraction for the L°° norm. 
We have the 


Lemma 5.1. Let be a £ R, s > 1 and p,pi , P 2 , q £ [l,oo] such that 

1 s 1 1 

- + - >- 1 -. 

P 3 pi p 2 


Then there exists £q > 0 and C > 0 such that for all u £ Bc(0,£o) and p £ 
Bp tQ ( R 3 ,R 8 ) nL°°(R 3 ,R 8 ), such that 

(5.1) 


\\(Q) a N(u, p)\\ Bs <c(\u\ + \\p\\ Ba ) \\v\\ L oo \\{Q) a v\\ B ° 

P.9 \ P2 >9 / Pi,9 

Proof. We recall the definition 

N(u, p) = VF(S(u) + rj)~ VF(S{u)) - d 2 F{S{u))p. 


We have 

N(u, p) = f f d 3 F(S(u) + d'ep)-p-9pdd'dd. 

Jo Jo 

Since for s £ R+, p, pi, P 2 , £ [l,oo] such that i + f > we have 

llitdls* < CIIitll||u||_B s ■ Then since s > 1 , we use (see 1EV971 Proposi- 

P.g Pi,q' P2 > 9 / 

tion 2.1]) 

\\d 3 Fm B s / q <C(s,F, Halloo) 

Then using Lemma rm we conclude the proof. 


Hence we have the 

Lemma 5.2. There exists 5o > 0 and £o > 0 such that for any 6 £ (0, <5o) 
and £ £ (0, £o)> f or an y u £ U{£,8) and z £ Z(u,S), Q u ,z(ot) maps 12(5) into 
itself. 

Proof. We have by means of Estimate o with e.g. a < —3, s = 0, p = q = 2 
and pi = P 2 = 4, if uq £ C and zq £ TL c {uq) D ff® (R 3 , K 8 ) are small enough 

OO | 

+ | {dS(u(v))A(u(v),p(v))(N(u(v), p(v)), dS(u{v))} JSf (u(v))) | 

+ \((dSj t (u(v)))A(u(v),p(v))(N(u(v),p(v)),dS(u(v))), JSf (u(v)))a ± (u)| 

+ | ((dSf {u(v)))A(u(v), p(v)){N(u(v),p(v)), dS{u{v))), JS^(u(v)))a T (v)\} dv 
< C6 2 {t)~ 2 . 

Hence for small S and small e, we have £7 U)Z (12(5)) C 12(5). 

To prove that Q u>z is a contraction for the L°° norm, we use the 
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Lemma 5.3. Let be a £ R, s > 0 and p,q £ [l,oo] such that sp > 3. Then 
for any e > 0 and M > 0 there exists C > 0 such that for all u, v! £ i?c(0,e) 
and 77, rj £ B^ q (M. 3 , R 8 ) n L°°(]R 3 , B r s(0, M)), such that 


11(0)" {N(u,p) N{u',r/)}\\ B ." < c| (\\(Qrv\\B^ q + \\(Q) ai v'\\B^ q ) 2 x 

x (lu-u'l + IKQ)^ {v-rf)\\ B . q ) 


+ l«l + |ti , | + 


(QV'V + (Q)^rf 

B s 

p ,q 

x (||(Q) </2 » ? I|b» >8 + \\(Q) a ' 2 v'\\B° p ^ ||(0)°« (v-v')WB- 

with 2o\ + a 2 = cr[ + ct' 2 + a' 3 = a. 

Proof. Since, we have 

N(u, v)= [ [ d 3 F(S(u) + O'dp) -p-dp dd'dd. 

Jo Jo 

we can also restrict the study to d 3 F((f) — d 3 F{(f'). If d 5 F ^ 0, we have 


\\(Qy{d 3 F(4,)-d 3 F(<j>’)) \\ B . iq 

< f nd 4 F(0+w- m Bi , \\(Qn<t >- n\B- q * 

Jo 

then since s > 0, we use 

\\d 4 FW)\\ B s q <C(s,F,M B ,J. 

Then using Lemma 14.11 we conclude the proof when d 5 F ^ 0. Otherwise the 
proof is easily adaptable since d 4 F is a constant matrix of Ad^C). 

We also need the 

Lemma 5.4. Let be a £ R, s > 0 and p,q £ [1, oo]. For any e > 0 and M > 0, 
there exists C > 0 such that for all u, u' £ -Bc(0,e) and p , p' £ Bp q (S. 3 , R 8 ) D 
L 00 (R 3 , -Brs(0, M)), one has 

\A(u,p) -A(u',p')\ < c{\u-u'\ + \\{Q) a {p - p}\\ B s^ 

Proof. We recall that 

A(u,p) = [{JdS(u),dS(u)) — {Jp,d 2 5(zt))] _1 . 

We have 


-A(zt, ?y) — A(u', p') = — [( JdS{u ), dS(u )) — (Jp, d 2 S(u)}] 1 x 
x {(JdS(u),dS(u)) - ( Jp 7 d 2 S{u )) - {JdS{u'),dS{u')) + (Jif,d 2 S(u'))} x 

x l(JdS(u'),dS(u')) - ( Jp l ,d 2 S(u'))}~ 1 . 
The lemma then follows from Lemma rm 
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Hence we have the 

Lemma 5.5. There exists So > 0 and e o such that there exists n £ (0,1) such 
that for any S £ (0, <5o) and e £ (0,£o) , for any u , u' £ U(e) and z £ Z(u,S) 
and z' £ Z(u\S), for any a, a' £ C(S), one has 

II Gu',z' ) — Qu,z{ol) || £oo( R +) 

< K (\u — 'u| Zj00 p R+ ) + |a — a| ioo ^ R+ ) + ||2: — z \\(R+,B^ 2 )) ■ 

Proof. It is a straightforward computation based on Lemma |5.3I with e.g. a < 
—6, cr2, <73 < —3 and s = 0, p = q = 2, on Lemma 15.41 on Lemma 15.11 with 
e.g. <7 < —3, p = q = 2 and pi = P 2 = 4 and on Lemma f4.31 

We now state the 

Lemma 5.6. There exists <5o > 0 and Eq > 0 such that for 8 £ (0, (5o) and 
e £ (0, £0) and any u GU(e) and z £ Z(u,8), the equation 

= E ± {u)a ± + ( JN(u , 77), JS±(u)) 

+ {dS(u)A(u, r]){N(u, 77), dS(u)) JSf 1 (u)) 

- {{dSf(u))A(u,r])(N(u,r]),dS(u)), JSfiu))^ 

- ((dSf (u))A(u, rj)(N(u, if), dS(u)}, JSf(u))a T ,, 

where 77(f) = a + (t)S+ (u) +a~(t)Sf(u) + z(t), has a unique solution in £2(5). 
Proof. The proof is now classical since we proved that the integral equation 

ait) = &u,z(a)(t) 

can be solved by means of the fixed point theorem. 


5.2. Step 2: Construction of z. Let be u £ U(e,8) and Zo £ H c (u(0)) D fJ* . Let 
us write Uao = lim u(t), we define T UtZo (z) by 


%, Z0 (z)(t) = e JtH ^z 0 

- [ e J ^ H ^P c (u(v))J{E(S(u(v)))-E(S(u O0 ))}z(v)dv 

J 0 

+ [ {d 2 F{S{u{v))) ~ d 2 F(S(u 00 ))} z(v)dv 

Jo 

+ f e J{t - v)HM P c (u(v))JN{u{v),p{v))dv 

Jo 

+ f e J< ' t ~ v ' ,H( ' u ^P c (u(v))dS(u(s))A(u(v), 77 ){N(u(v), 77(7;)), dS(u(v))) dv 
Jo 

— [ e J( ' t ~ v ^ H( ' u °°\dP c (u(v)))A(u{v), r](v))(N(u(v),r](v)),dS(u(v)))ri(v) dv. 


We have 
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Lemma 5.7. There exists <5o > 0 and £q > 0 and C > 0 such that for any S £ 
(0, (So), for any e £ (0,£o) and for any u £ U(e, 8), the application T u , zo maps 
Z(u, S ) into itself if \\zo\\ H s' < C6 

Proof. With Lemma rm and Lemma, T-l.41 we obtain 

WII H s ' 

^ < ^'\\^ > c{v, oo )T U}Z0 (t)\\ H s' 

< C\\z 0 \\ H s' +C f ||{E(S(u(v))) - £(S( Uoo ))} z(v)\\ Ha . dv 
Jo 

+C I ||{d 2 F(5( U (t;)))-d 2 ^(^ 00 ))}2(t;)|| ffs , dv 


+C [ \\N(u{v),r]{v))\\ H ^ dv 
Jo 


+C / \\dS(u{v))A(u(v),r](v))(N(u(v),ri{v)),dS(u(v)))\\ H .> dv 


+C / \\(dP c (u(v)))A(u(v),ri(v))(N(u(v),r](v)), dS(u(v)))r)(v)\\ HS ' dv. 

Jo 

Now, with Lemma ED we obtain 

\\%i,z 0 (,t)\\ H s' < C\\z 0 \\ H s' +Ce |u(v) - Wool IMI/fV dv 

Jo 

+C f (|u(t>)| + ||jy(v)|| H .') IkMILoc ||r?(t>)lljr«' dv 

Jo 

+C f (|u(t>)| + ||r?(v)|| H „') ||»?(w)||i«. ||»7 («)|Ih»' dv, 

Jo 

and so 

ll^“,2o(^)llff s ' — C|ko||_f/ s ' + Ce8 3 + C (e + (5) (5 2 + C [e + (5)” S 2 . 

Then, we also have 

r U) 20 (i) = e- ltH+1 ti E M r » dr zo 

rt 

+ / e -*t-v)H+iJ v EMr))dr pc ( u (v))Jd 2 F(S(u(v)))z(v)dv 

Jo 

rt 

+ / e ~ 1 ^ t ~ v l H + 1 fv S (“M) dr P c (u{v))JN(u{v), rj(v)) dv 

Jo 

+ [ e - I (t-v)J?+ I /*15 ( « ( r))drp c ( u ^))x 

Jo 

x dS(u(v))A(u(v),ri(v))(N(u(v), r?(t>)), dS(u(v))) dv 

_ f E(u(r)) dr x 

Jo 

x (dP c (u(v)))A(u(v), rj(v))(N(u(v),ri(v)),dS(u(v)))r](v) dv. 


52 


Nabile Boussaid 


Hence by Lemma ^31 and Theorem I I .21 we have 

ll”^u,Zo WII 2 

< C(t)~ 3/2 \\z 0 \\ B i3+3 + C [ (t - v)~ 3/2 \\d 2 F(S(u(v)))z(v)\\ B p+3 dv 

1,2 J o 1,2 

+C [ {t-v}~ 3/2 ||A^(m(w),77(w))|| b/3+ 3 dv 
Jo 1,2 

+C f (t — v)~ 3 ^ 2 \\dS(u(v))A(u(v),ri(v))(N(u(v),ri(v)),dS(u(v)))\\ b p+ 3 dv 
Jo 1,2 

+C f (t- v)~ 3/2 x 

JO 

x ||(dP c (u(u)))A(u(i;),r?(i;))(Ar(u(i;),r7(i;)),d5(u(v)))?7(v)|| s 3+3 dv. 

With Lemma o we infer 

IIW*)II B £, 2 — 

C(t) _3/2 ||2;o|| B 3+3+C [ {t - v)~ 3/2 \u(v)\ 2 \\z{v)\\ h / 3+3 dv 
1 ' 2 Jo 

+c [ {t- v)~ 3/2 (|u(u)| + 11^)11^3+3) lk(v)IL- h{v)\\ H0+ 3 dv 

JO 

+C [ (t- v)~ 3/2 (|u(u)| + ) ||77(u)|| i0 o h(v)\\ Hf3+3 dv 

Jo 

+C f {t~v)~ 3/2 {\u(v)\ + \\r]{v)\\ Hl3+3 )\\r)(v)\\ Lao ||t 7 (t ;)||^ + 3 dv. 

Jo 

With the estimate 


f (t- v) 3/2 {v) 3/2 dv<C(t) 3/2 , 

Jo 


we infer 

(t) 3/2 ||T u , Zo (7)|| B 3 2 < C\\zo\\ B ^ + Ce 2 6 + C(e + 6)5 2 +C{e + S) 2 S 2 . 
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Then we also have 


rt 

T u Zo (t) = e~ ltH+1 f °^ B (“( r )) dr zo + / e- 1 ( t ~ v ' >H+1 tt dr x 

Jo 


o 

x P c (u(v))JVF(r](v)) dv 

f e -i(*-o)H+i/*E(«(r)) dr x 


JO 

x P c (u(v)) J{X7F(S(u(v)) + rj(v)) — VF(S(u(v)) — VF(r](v))} dv 

rt 


0 —i(t—v)H+i f * E(u(r )) dr 


X 


/ 0 


x P c (u(i;))d5(u(u))A(u(i;), r](v))(N(u(v), r](v)), dS(u(v))) dv 

_ f e -i(t~v)H+if* E(u(r)) dr x 


x (dP c (w(v)))A(it(i;), rj(v))(N(u(v ), ri(v)), dS(u(v)))rj(v) dv. 


We now use Lemma E^l and Theorem ll.il except for the second term of the right 
hand side for which we used Theorem ll.2l since a > 3/2. We also use Lemma lo. II 
except for the third term of the right hand side for which an obvious adaptation 
of the proof of Lemma IQl gives 


II VF(S(u(v)) + V (v)) - VF(S(u(v)) - VF(t}(v))\\ h . 

< C(|u(u)| + ||t/(u)|| jjs ) |u(u)| ||r?(u)|| jF/3 


and so we obtain 


(t} 3/2 ||T u , Z0 (t)||^ 

< CUzollff. +CS 3 + C{£ + S)S + C(e + S) S 2 + C {s + S) 2 6 2 . 

Therefore we have that T U:Zo leaves Z(u, S ) invariant if ||^o||ij»' j <S and £ are small 
enough. 

An important property of T is given by the 

Lemma 5.8. There exists (>o > 0 and £o > 0 such that there exists k £ (0,1) 
such that for any 5 £ (0, (So), for any e £ (0, £o), for any u, v! £ h/(e,d), for 
any zq £ Tt c (u(0)), for any z r 0 £ Tt c {u'(0)), for z £ Z(u,6) and for any z' £ 
Z(v !, (S), one has 

— \\ z 0 ~ 2 ollz,~(R+,£f' ! ') + K |||w — W || L oo + \\z — Z |lioo( R + iffs ')| ■ 
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Proof. It is an easy consequences of straightforward estimates on the following 
identity 


'I'u.Zq (f) - ^ 


_ p-ltHiUoe) 


zo 


- [ cluivfiJiEiSiuivM-EiSiuooMzWdv 

Jo 

+ [ e- l{t - v)H(u ^P c (u{v))J (d 2 F(S(u(v))) - d 2 F{S( Uoo ))) z(v)dv 

Jo 

+ ( e - 1 (*-' l ’)ff(«oo)p c (, lt (, l ,)) jn(u(v), rj{v)) dv 

Jo 

+ ( e~ l( - t ~ v ^ H( ' y " X} ^P c (u(v))dS(u(v))A(u(v), rj(v))(N(u(v), r](v)), dS(u(v))} dv 

Jo 

— ( (dP c (u(v)))A(u(v), rj(v))(N(u(v), r](v)), dS(u(v)))rj(v) dv, 

Jo 

based only on Lemma ED 031 and ESI and on the fact that 
Pc(woo) P c (u(v)) = 

- E(S( Uoo )) e~ ls ' H ^ ds'P c (u(v)) 

+ Pc(« 00 )||e-’(- s ') ff (’ l “)(d 2 F(S(^)) 

- d 2 F{S{u 00 )))e-' s ' H ^A ds'P c (u(v)) 


P, 


form a family of operator in B(H S (R 3 ,C 8 ),H S (R 3 ,C 8 )) (to this end we use 
Lemma E3 such that 


Pc(«oo) (e~ iaH ^ - e~ lsH P c (u(v)) 


< Ce lu'oo - Wool , 


with C independent of u,u'. 


Lemma 5.9. There exists 5 > 0 and e > 0 such that for any u £ U(S,e) there 
is a solution z £ Z{u, 5) of the equation 


{ d t z = JF[{il)z + P c {u)JN(u, rf) 

+P c (u)dS(u)A(u, r])(N(u, rf), dS(u)) 
-{DP c (u))A(u, ri){N(u, rf), dS(u))r], 

z ( 0 ) = zq, 


where git) = a + (t)Sf (u) + a ( t)S 1 ( u ) + z(t), whenever zq £ -ff* is small 
enough. 


Proof. It is a consequence of the fix point theorem applied to T u ,z 0 
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Lemma 5.10. Under the assumptions of Lemma E3 for any u £ U (<5, e) and 
solution z of E3, with zq £ Hp small, the following limit 

Zoo = lim e ltH -‘ti E (ulr))dr z ^ 
t —►oo 

exists in H s (~l B ^ 2 D Ht a . Moreover, we have Zoo £ 7d c (0) and 

|| e -itJ?+i/ 0 ' E(u(r)) dr Zoo _ < cs 2 , 

\\e-' tH +'ti dr Zoo - z(t)\\ B , < C-f^, 

°°’ 2 w 

|| e - lt H + ^E M r)) dr^ _ 

Proof. Using exactly the same method as the one of Lemma rm applied to 


et* H - 1 ti E W) dr z(t) = z 0 

+ f e' vH ~'fo E «r)) drp c ( u ( w )) JJV(u(u), Jj(«)) dv 

Jo 


+ / e lvH ~ 1 Jo E< - u ^ dr x 

Jo 

x P c (u(v))dS(u(v))A(u(v), t](v))(N(u(v), r](v)), dS(u(v))) dv 

e ivH-iff E(u(r )) dr x 

x (dP c (u(v)))A(u(v), r}(v))(N(u(v),rj(v)), dS(u(v)))rj(v ) dv, 

we prove that the limit exist by the same way we also obtain the convergence 
rate. Since e~ ltH z x tends to zero, Zoo necessarily belongs to H c ( 0). 

Remark 5.1. The preceding proof also work with the formula 

pt 

e itZ? m -i/ 0 ‘ E(u{r)) dr z ^ = ~ Q+ e ™Dm.-iff E(u(r)) *-p c ( u ( v )) Vz ( v ) dv 

Jo 

+ [ e lvDm ~ l JJ E ^ u{r)) dr P c {u{v))J N {u{v),r]{v))d V 
Jo 

_|_ f e wD m -i ffE(u{r))dr x 

Jo 

x P c (u(v))dS(u(v))A(u(v), r](v))(N(u(v), ri{v)), dS(u(v))) dv 

e ivD m -if 0 ” E(u(r)) dr x 

x (. dP c (u(v)))A(u(v ), rj(v))(N(u(v),ri(v)), dS(u(v)))rj(v) dv. 
Hence we obtain the same result with 

e -itD m +if* E(u(r)) dr 
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instead of 

-itH+if* E(u(r)) dr 

c- ^oo • 

But we obtain the estimates 

E(u(r)) dr-^ _ ^ ||^ ( < ^ 

|| e - lt D m+ 1 f*E(u(r» dr Ioo _ Z (t)\\ B ^ < C-^, 

|| e -UD m+ 1 f 0 *E(u(r)) dr^ _ z ( t )|| ff ^ < C JL. 

5.3. Step 3: Construction of u. Here we want to solve the equation for u. We 
notice that z and a have been built in the previous section and are functions 
of u and z 0 £ H c (u(Q)). Let us introduce for any a £ S?(S) and uq £ -Bc(0, e) 
the function on U(e, 5): 

fu 0 (u)(t) = u 0 - [ A(u(v),rj(v))(N(u(v),r](v)),dS(u(v)))dv, 

Jo 

where r](t) = a + (t)Si(u) + a~(t)Sf(u) + z(t). We have the 

Lemma 5.11. There exists <5o > 0 and £q > 0 such that for any <5 £ (0, <5o)> 
for any £ £ (0, £o), the function f Uo maps U{e,5) into itself if uq and zq £ 
U* H H c (uo) are small enough. 

Proof. By means of Lemma rm we obtain 

\fu o (u)(t)\<\u 0 \ + C [ \\N(u{v),r]{v))\\ H s < K| +C(e + 6)S 2 . 

Jo 

Hence for uq and 5 small f Uo (u)(t) £ Bc(0,e). Estimate 15. Ill also gives the 
existence of (/ Uo (w)) = lim f Uo (u)(t ) and then 

t —>+00 

/ +oo /-I 

\\N(u(v),r](v))\\ H s_' r < -^(e + S)6 2 . 

The function f Uo has also a Lipshitz property as stated by the 

Lemma 5.12. There exists So > 0 and £q > 0 such that there exists k £ (0,1) 
such that for any S £ (0,<5o), for any e £ (0,£o), for any u, u’ £ U{e,5), for 
any zq £ Tt c (u{0)) D Lf* , for any z' 0 £ H c (u'(0)) D H * small enough, for uq , u' 0 
small enough, there exists k £ (0,1) such that 

\fno(u) - fu< Q (u') | ioo(R+) < | Wo - u'o\ + K (ll u - u 'IIl°°(R+) + \\ z 0 ~ z 'o\\Hf) • 
Proof. This a straightforward consequence of Lemma 15.3115.4115.51 and EH 
We are now able to prove the 
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Lemma 5.13. There exists So > 0 such that for any S G (0, <5o) such that for 
any uq € C small and zo G 7t c (zio) D Hf, small, the equation 

fit =-A(u,r]){N(u,T}),dS(u)), 

\ u{ 0 ) = u 0 , 

where 77 (f) = a + (t)Si (zt) + or (t)Sf (u) + z(t), has a unique solution in U(5, e). 

Proof. This is also a straightforward consequence of the fixed point theorem 
for fu 0 - 


5.4. Step 4 : End of the proof of Theorem roi We now conclude our proof with 
the 

Lemma 5.14 (Decomposition lemma 2). Let be s > 0 and p > 1 there 
exist 6 > 0 and a C°° map Uq ■ Bw s ’p(0, S) 1 —> l?c(0,£) which satisfies for if G 
B W s,p( 0 , 5 ) 

if = S(u) + Tj, with 77 G (</>o}‘ L •<=>■ u = Uo{if) 

Proof. In fact, we just write if with respect to the spectral decomposition of H: 

if = u(fo + r = S(u) + 77 

where r G {<^ 0 }^ and 77 = r — h(u), with h defined in Pr 0 D 0 siti 0 n ll.ll 

With respect to the notation of Theorem 11.31 for ifo = S(uq) + a + (0)Si(uo) + 
a _ (0)S'j _ (uo) + P c (uo)5o where z~o G Ran(P c ), we introduce 

fu 0 = U 0 (tfo) 

Uo = Pc (ifo ~ S(vo)) 


(G(u 0 ,z 0 )i = U 0 (S(u 0 ) + a + (0)S'j f (u o ) + a (0 )S 1 (zz 0 ) + Pc(wo)^o) 
J G(u 0 ,z 0 )2 = P c (s{u 0 ) - S{G(u 0 ,zb)i) 

[ +a + (0)S , j f (zto) + a _ (0)5'j _ (zi o ) + P c (zz 0 )-So) 


Then using Uq(S(uo)) = uq, we write G(uq , zo) = (uq, zq) + G(uo, zo ), with 
G(u 0 ,z 0 ) - G{u'o,zo) < k (|zt 0 - Uq| + ||zS - zo\\ HB ,^j 


with k < 1/2 if zto, u' 0 and Zq, zq small enough. Hence in this case G is invertible 
with a Lipshitz inverse F. Then we choose 


E{vo,t;o) 

= {S(F(v 0 ,£o)i) 

+a(F(v 0 ,fo)i) + (0)Sf(F(v 0 ,£ 0 ) 1 ) + a(F(v 0 ,fo)i)~(0)Sf(F(v 0 ,fo)i) 
+P c(F(vq, £o)i)F(v 0 , £0)2 - S(v 0 ), (/ 1 ) 
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and 

and 


Coo = (P c (F(^o,Co)l)^(^0,Co)2) c 


Eoo = / {E(F(v 0 ,Zo)i(v))-E((F(vo,Zo)i) 00 )} dv. 

Jo 

In the proof of Lemma rm we see that S is of the same order as ||Co||flV- The 
rest of the Theorem easily follows. 

□ 


APPENDICES 


A. The wave operator and similarity for the linearized operator 

Inspired by |Katfifij . we use an argument of similarity to prove the 
Lemma A.l. For all s £ R + , there exists C s > 0, such that 

Vf eK, ||e tJH ^\\ C {H‘)<C s . 

We prove this lemma by using the boundedness in H s of the wave operator: 
W± = s — lim e -tH( 2 rJ e -it(H-E(z))p c ( H } 

t —t±oo 

and the intertwining property: 

e -mzYJp c ( z y = w ± e- lt( ' H - E( ' z) '>P c (H)(W ± y g , 

This boundedness follows from the 

Lemma A.2 (Smooth and small non-selfadjoint perturbations). Let be 

ip £ L 2 and a > 1. Then there exists e > 0 and C > 0 such that 

POO 

\/z£B c (0,e), || <Q>-° e sJH ^P c (z)iP\\ 2 2 ds<C\\iP\\l (A.l) 

Jo 

Proof. By Lemma 14.41 . we prove P c (z)R{z, 0)P c (0)P c (z) = P c (z). Taking the 
adjoint with respect to the real structure, we infer P c (z)P c (0)i?(z)*P c (z) = 
P c (z). Then, we write 

||(Q)-V J *WP C (*)|| 

= C\\(Q}-°P c (z)e tJH ^P c (0)R(z)*P c (z)\\ 

< \\(Q)- a P c (z)e-^ H - E ^P c (0)R(zyP c (z)\\ 

+ [\\(Q)~ a P c (z)e l ' t ~ s ^ JH< ' z ^D'VF(S(z))e~ ls< ' H ~ E( ' z ^P c (0)R(z)P c (z)\\ ds 

Jo 

< \\(Q)-^P c (z)e-^ H - E ^P c (0)R(z)P c (z)\\ 

+C\z\ 2 [ ||(Q)- CT P c ( 2 )e( t - s ) J ^(Q)- CT ||||(O)- CT e- ls ^- B ( z »P c (0)|| ds. 

Jo 

Using Proposition Emu we obtain the claim ED for 2 sufficiently small. 
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This give us the existence and the boundedness of the wave operator, as stated 
by the following 

Lemma A. 3. Let be Wt = J e -tt(H-E(z))p C (H). Then the limits 

W ± = lim Wt 

t —»=hoo 

exist in B(H S ) and their range is Ran (P c (z)). The same is true for Wf and 

(IV±) -1 = lim (W t ) -1 . 

' ' t —► ioo 

Proof. Let us define Wt = e~ tH ^ zS> J e - lt ( H - E ( z ))_ we have for p £ H c (z) and 

e h c { 0) 

(cf, W t ip) = (</>, ip) + /</>, ^ W S V^ ds, 

Since we have 

< C|z| 2 || (Q) -17 e tJH( ' z ' ) cj)\\\\(Q)~ ,J e~ lt( ' H ~ E( ' z ^ip\\. 

which gives (p, j^W s tp) £ L 1 (M), so W± exists and is bounded in C{H c (f)),H c (z)) 
by the previous lemma. Since for any vector p in an eigenspace of JH(z), Wfcp 
tends weakly to zero, we obtain that the range of W ± is a subspace of the range 
of P c (z). Then the same statements about (Wt ) -1 follows by the same way. The 
invertibility is then immediate. 

Proof fof Lem,m,a, l A. 1\) . The L 2 bound follows from the intertwining property as 
explained before Lemma fA.21 

The proof of the H k bounds follows from commutation argument, we apply the 
same scheme to 

d t e- tJH ^P c {z) = [^,P c ( 2 ;)]c- tJff Wp c ( 2 !) + P c («)c- tJff W[^,P c {z)\ 
+P c (z)[d i ,e- tJH ^]P c (z) 

= [d u P c (z)]e- tJH ^d i + P c (z)e- tJH ^d u P c (*)] 

+ [ e- {t - s '> JH{z '>P c {z)(d i D\/F(S))e- sJH ^P c {z)dz. 

Jo 
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